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I . INTRODUCTION 


Thli grant covared a period of nore than twelve yeare and aupported a 
wide variety of ionoepheric particle and field research performed with sound- 
ing rockets. This research can be grouped into seven major rocket projects. 

Each project included the usual phaeee of experiment design and construction, 
flight operations, and data analysis. A diecussion of the experimental de- 
tails and scientific results of each project is given in section II. 

Several important experimental techniques have been developed uuder sup- 
port of this grant. The double-probe electric field technique was perfected 
during this series of experiments and the mechanical design of the electric 
field booms is now widely used by other experimenters as well. Other Instru- 
mentation developed under this grant includes swept frequency receivers, density 
fluctuation detectors, ion drift detectors and numerous charged particle de- 
tectors and iiiass spectrometers. 

Most of the work done under this grant has been carried out by a succes- 
sion of graduate students who used the results for Ph.D. theses. The students 
are M. C. Kelley, C. W. Carlson, E. A. Bering, R. B. Torbert, M. A. Temerin, 
and A. J. Mallinckrodt. 

II. TASKS PERFORMED AND SCIENTIFIC RESULTS 

A. Javelin NASA 8.55, 8.56 (Auroral Zone, 1970 ) 

Two Javelin rockets were launched from Fort Churchill, Canada during the 
spring of 1970 to investigate the electric field and particle precipitation 
structure of discrete auroral arcs. The instrumentation of these payloads in- 
cluded double probe E-field detectors, an ion drift detector, energetic elec- 
tron and ion electrostatic analyzers, a split Langmuir current probe and 
several photometers. Of these Instruments, the double probe, the ion drift 
detector and the split Langmuir probe provided three independent measures of 
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•lactric fields. The results verified thet thesv measuring techniques, which 
wure somewhat unproven at the time, gave consistent results (Moser et al . , 
1971T). (Note: References with "T" are listed in "Talks Presented at Scien- 

tific Mcetinp^.") The operation and results of the ion drift detector are 
given by Carlson (1971T) and Carlson and Kelley (1977). The scientific re- 
sults of these flights arc discussed by Keli.ey at al . (1975), Bering et al . 
(1975), Kelley and Carlson (1977), and Carlson and Kelley (1977); and in the 
Ph.D. thesis of C. W. Carl eon entitled "Rocket measurements of auroral zone 
low energy charged pat. rides" (University of California, 1974). Results from 
these flights also motivated the study of electric field and field aligned 
currents discussed by Mallinckrodc and Carlson (1978). 

B. Black Brandt III. NASA 2201. 2202 (Polar Cap. 1971 ) 

Two Black Brandt III rockets were flown from Resolute Bay N.W.T. in 1971 
to investigate the mapping of electric fields and charged particles from the 
earth's deep magnetotail down to the polar cap ionosphere. These experiments 
were made in conjunction with the lunar orbiting UC subsatellite released from 
Apollo 15. Results from this experiment are presented by McCoy ct al . (1976). 

C. Dual Hawk (Equatorial Spread F, 1972) 

The Dual Hawk rocket project was a collaboration between UC Berkeley and 
the Max-Planck-Institut fUr Extraterrestrische Physlk in Garchlng, West Ger- 
many. The objective of this experiment was to trace the development of equa- 
torial Spread F using a combination of optical observations of barium releases 
and In-sltu measurements of plasma density and electric fields. Our vole in 
the experiment was to provide the in-situ instrumentation package. These 
flights were unsuccessful owing to launch vehicle failures. 

D. Javelin NASA 8.62, 8.63 (Equatorial Spread F, 1973 ) 

Two Javelin rockets were launched from Natal, Brazil In 1973 to investi- 
gate the origin and evolution of equatorial Spread F. These flights were made 
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in collaboration with tho Max-Planck-Inatitut and included barium releaaea 
in addition to the in>aitu meaaurementa. Improved faat Langmuir probea auited 
to measure electrostatic waves and irregularities were developed for this ex> 
periment (Kelley at si . . 1975) • The scientific rci>u<.ts of these flights were 
evidence for a Rayleigh-Taylor instability during Spread F and the first direct 
evidence for plasma "bubbles" or buoyant low-density regions appearing during 
Spread F conditions (Kelley et al .. 1976). 

E. Polar Cusp Black Brandt IV (Greenland. 197-75 ) 

The Greenland polar cusp campaign was a collaboration between UC Berkeley, 
the Max-Planck*-Institut fC(r Exterrestrische Physik and the Danish Meteorolog- 
ical Instltut. We supplied instrumentation for two rocket payloads which in- 
cluded E-fleld double probes, a plasma density experiment and an array of ion 
and electron detectors. The payloads also contained shaped-charge barium ex- 
periments supplied by the Max lanck-lnstltut. The main objectives of this ex- 
periment were to Investigate the plasma flow pattern near the cusp, study sola 
wind plasma entry into the cusp, and search for magnetic field aligned electric 
fields and particle acceleration. 

Both flights were very successful in fulfilling these goals. The flights 
took place inside the dayslde auroral region lying Just equatorward of the con- 
vection reversal boundary (Jorgensen et al . , 1980) . Although this region is 
nominally on closed magnetic shells, we found unambiguous evidence of solar 
wind plasma injected onto these field lines. These injections consist of brief 
bursts lasting a few seconds and separated by periods of about 100 seconds. 

This property of the injection allowed us to use time of flight analysis to 
ascertain the ion mass composition and verify that the injections were of 
solar wind composition (Torbert and Carlson, 1976; Carlson and Torbert, 1980). 
These flights also provided evidence for magnetic field aligned acceleration 
of particles. This evidence came from the barium shaped charge and from the 
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low energy electron detectors (Torbert end Cerlu.on, 1980). Intereetlng new 
tMeeurenente of low frequency waves were also obtained (Ti.aierin and Kelley, 

1980; Temerin and Parady, 1980). 

These flights were thesis projects for two graduate students, M. A. Tea- 
erin and R. B. Torbert 

F. Project Porcupine (Auroral Zone, 1977. 1979 ) 

Project Porcupine included four rocket flights over a period of four years. 
It was organized and managed by the Max-Planck-Institut and included several 
other European groups in addition to the U.S. collaboration from UC Berkeley. 

The payloads consisted of a main "mother" payload which included extensive 
particle and field instrumentation and four ejected subpayloads. For each 
payload, wc supplied two subpayloads which contained R-fleld probes, Langmuir 
probes, a magnetometer, and a complement of particle detectors. The other two 
cannisters contained a barium shaped charge and a xenon plasma gun. These sub- 
payloads were ejected from the main payload during flight and allowed three- 
point measurements with separation distances up to four kilometers. The launch 
vehicles were "Arles" rockets. The first launch had a vehicle fa.'lure, but the 
remaining three flights were successful and have returned valuable data. 

Preliminary results from these flights have been presented in several sci- 
entific meetings (Carlson et a l. , 1977T; Malllnckrodt and Carlson, 1977T; Mal- 
llnckrodt et al ., 1978T; Carlson, 1979T; Ergun et al . , 1979T; Malllnckrodt, 
1979T; 'sudeki et al . , 1980T; Witt et al . , 1979T; Bush et al . , 1979T; and Bush 
et al . , 1980T) . One of the unique observations on these flights was strong 
emiF.slon of Bernstein waves near the upper hybrid frequency. These waves 
reached amplitudes of tens of mV/m. In one case the wave amplitude was suf- 
ficient to produce parametric decay into other lower frequency modes. The 
barium release also produced unexpected effects upon the wave propagation, 
causing trapping of waves inside the enhanced plasma of the barium cloud. The 
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other unique opportunity afforded by theao fllghte wae to study effects of the 
xenon plasna gun. Thie analysle is still under way and preliminary results 
were presented by Bush et al . (1980T). First results of these flights are en- 
closed In the Appendix and the first three publications are In preparation. 

One Ph.D. thesis, by A. J. Malllnckrodt, was based upon these experiments, and 
a second, by R. Bush, Is In progress. 

G. Project Flrewhael (Active Magnetospherlc Satellite. 1980) 

Project Flrewheel was an attempt Inject an artificial plasms cloud Into 
the tenuous magnetotali plasma and study Its effects such as magnetic cavity 
formation, magnetic diffusion, momentum exchange, etc. These objectives are 
discussed in great detail in the experiment proposal. 

The UC participation In Flrewheel was to supply a completely Instrumented 
subsatelllte to be ejected In orbit and perform In-sltu measurements of tha 
natural and artificially produced plasmas. The spacecraft was launched on the 
second Arlane test vehicle, which unfortunately experienced a first stage motor 


failure. 
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ROCKET OBSERVATIONS OF NATURAL AND 
ARTIFICIALLY STIMULATED PHENOMENA WITHIN A 
QUIET AURORAL ARC 

Albert John MaUinckrodt III 


ABSTRACT 

An Aries sounding rocket launched from Kiruna, Sweden measured particle fluxes, DC 
and AC electric fleids, and plasma density within and on the equatorward side of a quiet auroral 
arc. A striking anticorrelation between the cross>arc component of the electric field and the 
energy of the secondary peak in tl:e election differential particle flux was observed. A numeri- 
cal procedure has been developed to simulate auroral electrodynamics and shows the observed 
anticorrel.v'jn to be consistent with a model in which field-aligned sheet currents are negligible 
with respect to ionospheric height integrated currents. The electric field is then determined 
from the divergenceless flow of the ionospheric current system and the transverse conductivi- 
ties which arc in turn controlled by the precipitating electrons. In addition large amplitude high 
frequency plasma waves were observed within a plasma density enhancement artificially created 
by u barium shaped charge experiment included Obi the payload. The waves are associated with 
the formation of a plateau in the formerly peaked auroral electron distribution function. These 
observations arc explained in terms of the refractive trapping of wave packets within the 
artificial density enhancement. The resulting prolonged interaction of these waves with the 
auroral electrons allows the growth of the waves at the expense of the free energy in the elec- 


tron distribution function. 


Table of Contents 


ACKNOWLEDGEMENTS i 

Table of Contents ii 

I. INTRODUCTION I 

II. INSTRUMENTATION 8 

■ ■ A. Panicle enperimenis .10 

B. Electric ndit ex|tcriment!S 14 

C. Plasma experiments 15 

D. Magnetic field experiment 17 

III. LAUNCHES 18 

IV. DATA 21 

A Overview 21 

B. Anticorrclutiun of field and precipitation 2S 

C. Electron distributio'' function 25 

D High frequency plasma waves 27 

E. Summary 28 

V. ANALYSIS 29 


TtMn *f I'aMfM* 


iii 





A. On itw control of the tonovpheric electric Held... 29 

D. On the high amplitude wave oltservationi and the wavc*particle effects 43 

VI. CONCLUSIONS 58 

VII. APPENDIX • Numerical simulation of auroral electrodynamics 60 

A. Assumptions r the routine ....60 

B. Data flic structure 61 

C. Description of the routines 62 

f 

D. Examples 68 

REFERENCES 71 

TABLES 73 

FIGURE CAPTIONS 76 

FIGURES 78 



I 


I. INTRODUCTION 

During the past two decades the space research effort has revealed the large-scale struc- 
ture of that portion of space which is controlled by the earth’s magnetic field - the moKnefo- 
sphere. Prior to the first data returned to earth from artificial satellites, the region of space sur- 
rounding the earth was generally envisioned as a listless, cold, and extremely tenuous medium. 
The discovery of the radiation belts (Van Allen et ai, I9S8) provided dramatic evidence to the 
contrary and resulted in an awakened scientific interest in near-earth space phenomena. With 
the discovery of the solar wind, a continuous stream of plasma flowing from the sun, an energy 
source for the observed high energy particles was found and it was realized that the description 
uf magnetosphcric processes involved the interaction of a fully ionized solar wind plasma with 
the magnetic field of the earth. 

Until the early I970's the bulk of experimental work was directed toward gathering a body 
of data from which a comprehensive empirical model of the earth's magnetosphere could be 
established. While there is still some disagreement regarding the interpretation of the data, a 
reasonable version of the current magnetosphcric model is displayed in Figure I.l. The mag- 
netic field of the earth is confined to a comet-like region of space about 40 /(,. in diameter and 
perhaps 1000 /{,. in length. The blunt nose of the magnetosphere is formed by the dynamic 
pressure of the solar wind which compresses the magnetic field on the day side of the earth and 
confines its flanks further downstream. At the furthest downstream points the magnetosphere 
presumably becomes virtually indistinguishable from the interplanetary medium. To date only 
a few encounters with the tail have occurred beyond lunar distances. The boundary of the 
earth's m ignetospherc is called the magnetopause and is characterized by the transition from a 
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region of large plasma pressure outside the magnetosphere to a region of large magnetic field 
strength inside. The discontinuity in the magnetic field requires large sheet currents to flow on 
the magnetopause. 

In discussing the magnetosphere it is helpful to define three distinctly different typ<;;i of 
magnetic field lines. Closed magnetic field lines are those which pass only through the earth. 
Open magnetic field lines are those which pass both through the earth and the sun. Finally, 
interplanetary field lines are those which do not pass through the earth. The magnetosphere 
proper consists of closed field lines and the near earth portions of open field lines. Just outside 
the magnetosphere lies a region called the magnetosheatli. The magnetosheath consists of solar 
wind plasma on interplanetary field lines, but which differs from the undisturbed solar wind by 
virtue of its hav.ng passed through the how shock, the shock wave formed by the pile up of 
solar wind plasma at the nose of the magnetosphere. In this region the plasma is hotter and 
more dense and the flow is more disordered than in the solar wind particularly at the sub-solar 
stagnation point on the magnetopause. As we continue tailward the distinction becomes less 
clear. 

Satellites instrumented with particle and field detectors have mapped out the various 
regions within the magnetosphere and delineated the borders of these regions. Due to the adia- 
batic motion of the plasma these boundaries can often be defined in terms of the magnetic lati- 
tude of the field lines which are coincident with them. In all of the following discussion it is 
helpful to keep in mind that whenever the component of the electric field parallel to the mag- 
netic field vanishes the frozen-in condition applies. That is, one may think of the field and the 
plasma moving as a single entity. Due to the collisionless nature of the plasma throughout 
most of the magnetosphere, the conductivity along magnetic field lines is very high and as a 
result tlte frozen-in condition applies. In certain regions however, wave turbulence may result 
in an anomalous resistivity which can produce magnetic field-aligned electric fields and invalidate 
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At the lowest latitudes and thus close to the earth, the magnetic held lines lie close to 
their expected dipolar configuration. This region is called the plasmasphere. The plasma here 
tends to be stably trapped for long periods and diffusion of ionospheric plasma along magnetic 
field lines fills the plasmasphere to high densities. The phenomenon of shell-splitting due to 
the longitudinally asymmetric distortion of the magnetic field lines complicates this picture and 
results in the formation of a quasi-trapped region in which certain particles, depending on their 
magnetic latitude and mirror altitude are not able to complete an entire orbit around the earth 
(Stone, 1963; Anderson ei al„ 1965; Roederer, 1967; Roederer, 1970). The plasma here coro- 
tates with the earth due to the presence of an electric field generated in the ionosphere by the 
viscous coupling of the ionospheric plasma to the neutral atmosphere. Due to the high conduc- 
tivity for cut rents along magnetic field lines, the field lines may be considered equipotentials. 
Thus the ionospheric electric field is mapped out into the magnetosphere. Were it not for the 
presence of the solar wind or other outside influence's the plasmasphere would extend to arbi- 
trary distances and the oniy limit on its corotation would be the limited angular momentum of 
the earth available to be imparted to the plasma. The solar wind however imposes its own, 
quite different requirements on the motion of the plasma at the boundaries of the magneto- 
sphere and it is this competition which results in the formation of the other regions of the mag- 
netosphere. 

At a magnetic latitude of about 65 degrees the corotation stops (Mozer et aL, 1980). This 
cofotation boundary separates the plasmasphere from the plasma sheet. In the piasma sheet the 
field lines are still closed but highly distorted on the night side into a long tail with fieid lines 
running sunward to the north of and anti-sunward to the south of a neutral sheet dividing the 
niagnetotail into two regions of oppositely directed magnetic field. A neutral sheet current flows 
from dawn to dusk in the neutral sheet as can be seen from the tangential discontinuity in the 
magnetic field at this boundary. At the inner edge .'f the neutral sheet the ring current, 'm occa- 
sionally intense current, flows. Convection in the plasma sheet is predominantly sunward 
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owing to the proce» of magnetic held line reconnection which depletes the magnet, c flux on 
the day side and builds it up on the night side. This sunward convection continues around the 
dawn and dusk sides of the piasmasphere to the nose of the magnetosphere where it reconnects 
with the interplanetary held. The competition between the sunward and anti*sunward flow at 
the corotation boundary on the dusk side yields a simple determination of the corotation boun- 
dary position through satellite electric held measurements. On the dawn side the distinction is 
not so ciear since the flows are in the same direction. 

The plasma sheet field lines reconnect to interplanetary fl!*!d lines at the nose in the entry 
layer. On the day side a magnetic cusp is formed which allows direct penetration of solar wind 
plasma to the ionosphere. Further tailward the solar wind plasma, having mirrored in the cusp, 
is swept back on these now open field lines forming the plasma mantle just within the polar 
magnetopause. Due to the combined effects of equatorward convection and velocity dispersion 
of the tailward flowing mantle plasma, the equatorward portions of the mantle consist of ihe 
lowest energy particles (Paschman, 1979). The field lines convect into the polar lobes of the 
magnetosphere where the primary source of plasma is the polar wind due to diffusion of plasma 
from the polar ionosphere. On a given field line the boundary between the plasma mantle and 
the polar lobe is indistinct but it will likely progress tailward as we consider field lines lying 
closer to the equatorial plane. 

Along the flanks of the magnetosphere lies the boundary layer. In this region solar wind 
plasma is injected inside the magnetopause and diffuses onto closed field lines creating a viscous 
transfer of momentum and energy to the magnetosphere and causing anti-sunward convection 
(Lemaire, 1977; Lemaire and Roth, 1978). On the inside edge of the boundary layer, a similar 
viscous interaction with the plasma sheet plasma causes a sunward convection closing the circu- 
lation pattern. 

The picture described in the preceding paragraphs applies to highly idealized, steady-state 
conditions. Since the primary energy source driving magnetospheric phenomena is the solar 



wind, the validity of 'he model depends predominantly on its parameters. The orien.ation of 
the interplanetary magnetic held, the solar wind plasma density, temperature, and bulk flow 
speed, and the temporal variability and history of all these affect the state of the magneto- 
sphere. 

Before the advent of electric field experiments capable of operating in the tenuous plasma 
environment of the outer magnetosphere, researchers relied on plasma density measurements 
to define the boundary of the plasmasphere the phsmapause. The plasma density in the 
plasma sheet is reduced by about two orders of magnitude from that within the plasmasphere 
due to the loss of plasma through reconnection at the nose and the much larger volume of 
space occupied by plasma sheet field lines. In periods of long solar wind stability this method 
yields satisfactory results, but during an increase in the dynamic pressure of the solar wind, the 
magnetosphere is compressed on the day side and the corotation boundary is driven to lower 
magnetic latitudes. Plasma at the outer boundaries of the plasmasphere, built up by long 
periods of undisturbed diffusion of ionospheric plasma, is lost. When the corotation boundary 
moves outward, the field lines refill with ionospheric plasma on a time scale of days. During 
periods of high variability the plasma density at latitudes near the corotation boundary is a com- 
plicated function of the actual position of the corotation boundary over the last few days 
because of the long time scale for the ionospheric diffusion. 

Another effect of quite some importance is the orientation of the interplanetary magnetic 
field. Depending on this orientation, magnetic field reconnection at the nose and thus plasma 
sheet convection may be enhanced or inhibited. 

Rapid changes in any of the solar wind parameters induce changes in the state of the mag- 
netosphere and often result in the triggering of a masnetospheric subsiorm. The substorm is 
characterized by the energization of plasma sheet particles at the expense of the stored magnetic 
energy of the magnetotail. The high latitude portion of the plasma sheet includes the auroral 
:om>, those field lines on which electrons are energized and driven into the ionosphere creating 
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the most dramatic visual manifestation of the solar wind*magnetosphere nteraction • the 
aurora. 

The subject of this thesis is conflned to those processes occurring on auroral held lines 
and, more specifically, to observations of natural and artificial phenomena occurring within the 
ionosphere in the presence of auroral electron precipitation. For many years prior to the advent 
of artiflcial satellites, instrumented balloons had monitored high altitude phenomena in the 
auroral zone. The difficulty with the interpretation of these data was that the measurements at 
balloon altitudes represented averages of the ionospheric phenomena over large distances. For 
example electric held measuremei^ts at balloon altitudes may be considered to represent an 
average of the E*region field over a distance of about 100 km. The visual thickness of auroral 
forms is on the order of kilometers (Maggs and Davis, 1968) and thus much of the expected 
fine structure is lost. Furthermore the actual particles producing the displays and their energies 
could only be guessed at from spectroscopic observations of auroral emissions and models of 
the high altitude atmosphere. In addition, with only the data from ground based magnetome- 
ters to work with, the three dimensional current system associated with auroral displays was not 
well defined. 

Satellites have yielded valuable information on large scale structure of the auroral zone 
and the distribution of field-aligned currents, but due to their high velocity they give essentially 
no information on the fine structure or temporal rariations in auroral cisplays. 

The sounding rocket is uniquely suited to the needs of auroral investigations due to 
several factors. They can be launched on short notice into auroral conditions of the scientist’s 
choice. They permit relatively long, continuous monitoring of the conditions within a single 
auroral arc. They have very low costs both in terms of manpower and money compared to 
satellites. 


The data on which this research is based are from the second flight of the Aries-Porcupine 
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campaigns and. more speciftcally, from an ejcctable sub-payload designed and built as a portion 
of this thesis research at the U.C. Space Science Laboratory. Two very interesting results were 
obtained from these data. 

First, a very strict anticorrelaiion between the D.C. electric field and the peak energy of 
the precipitating electrons was observed. This observation will be discussed in terms of the 

'.MiUiiUCA 

influence that the energetic electrons exert on the ionospheric conductivity due to the ioniza- 
tivin they produce. A numerical simulation of the ionospheric electrodynamics was developed 
and the results of its application to these data will be presented. 

Second, high amplitude plasma waves were observed for a period following apogee during 
which the payload was enveloped by an artificial plasma density enhancement produced by a 
barium shaped charge experiment. At the same time a dramatic change in the distribution 
function of the auroral electron beam was noticed. The distribution, which displayed a prom- 
inent peak at an energy of 2-3 keV prior to the barium release and the observation of the 
waves, was modified to a plateau type distribution afterwards. These observations will be dis- 
cussed in terms of the trapping of upper hybrid wave packets by the artificial plasma cloud. 
These trapped wave packets, due to their prolonged interaction with the auroral beam, attain 
substantially increased growth. As a result they grow to large magnitude and the increased 
interaction with the auroral electrons depletes the free energy represented by the peak in the 
distribution function. In addition a parametric decay process appears to have been driven by 
the high amplitude upper hybrid pump. Possible decay modes and wave vector matching 
schemes for this process will be discussed. 

in all, four Aries-Porcupine payloads were built and launched. The first, launched in 
March 1976, yielded no scientific data as a result of the failure of the launch vehicle. The 
second was successfully launched one year later. The final two were launched two years later in 
March 1979 and will be the subject of other reports. 
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II. INSTRUMENTATION 

The ArieS’Porcupine project was an international collaboration involving researchers from 
several European countries, the U.S., and the Soviet Union. The use of the Aries rocket with 
its very large payload capacity permit led the construction of one of the most comprehensive 
instrument packages ever launched on an auroral sounding rocket mission. The main payload 
included the following experiments: 

1 Double probe electric Held 

a D.C. channel < 80 Hz 
b A.C. channel SO Hz < 20 kHz 

2 Plasma density fluctuations 

3 Low energy electrons 

4 R.F. probes 

a High frequency .8 - I4.S MHz for electron density and temperature measurements 
b Middle frequency .1 -1.5 MHz for Held aligned electron current measurements 

5 Three-axis fluxgate magnetometer (resolution ~ I nT) 

6 Search coil magnetometer for m.t^vteiic fluctuations 10 Hz - SOO Hz 

7 Suprathermal electrons 100 eV -25 keV 

8 Suprathermal protons I - 25 keV 
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9 Energetic particles 

a Electrons I S - 200 keV 
b Ions 2S keV - 2 MeV 

In addition, four free-flyi'v* sub-paytoods were included on each flight. These were 
released early in the flight and were used both for active plasma experiments and as indepen- 
dent and spatially separate probes of local plasma conditions. On the four Porcupine payloads 
flown during the campaigns, various configurations of these sub- payloads were included. The 
second flight, which is the subject of this research, was equipped with a barium shaped charge 
experiment, a ce.sium ion gun experiment, and two ejectuble multiply instrumented probes. 

The buiium shaped charge experiment created a large artilkiai plasma cloud near apogee 
due to the photo-ionization of barium atoms vaporized in the detonation of a mixture of 
barium and high explosive. Optical transitions in the ionized and neutral components of the 
barium allow observation of the clouds from the grou id and their motions can yield important 
information about neutral winds and electric fields both transverse and parallel to the ambient 
magnetic field. In addition the auroral plasma perturbation cun ariiflcially stimulate plasma 
processes in the ionosphere which may then be monitored by the instruments of the main pay- 
load and the cjectablc probes. Since the barium cloud must be illuminated by the sun and the 
observation stations must be in moonless dark, the inclusion of a barium shaped charge experi- 
ment limits the launch window to a few hours after sunset or before sunri.se a couple of weeks 
a month. 

The cesium ion gun was designed to yield an intense curr '* low energy cesium ions. 
Once again the clTccts of this major pciiurbaiion on the auroral plasma was to be monitored by 
the instrumented payloads. Due to the apparent malfunction of this device it will receive no 
further discussion. 

The instrumentation aboard the ejectuble probes was designed and built at the University 
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of California Space Sciences Laboratory. This instrumeniaiion provided a comprehensive set of 
plasma diagnostics. The experiments may be roughly classified os particle, electric field, 
plasma, or magnetic field. Table I lists the experiments along with their ranges of operation, 
sensitivities, and sample times. Figure II. I is a sketch of an ejectable probe along with detector 
aspects and payload coordinates. 

11. A . Particle experiments 

The payload was equipped with two different types of particle experiments which meas* 
c:ed low energy electrons and high'energy electrons and ions. 




Low energy electrons 

Retarding potential analyzers (RPA's) were used to monitor electron fluxes at energies 
from 6 eV to 500 cV. A diagram of this instrument is shown in Figure II. 2a. A detailed 
description of this detector is to be found elsewhere (Carlson, 1974). Briefly the instrument 
uses a repelling grid held at + 1 2 volts relative to the payload potential to eliminate the large 
thermal ion fluxes. Then a swept potential barrier is used to deny entrance to all electrons 
below a certain energy. Electrostatic focussing then allows only electrons of minima! excess 
energy (about 30 eV) access to the m.outh of a channel electron multiplier (CEM). The CEM 
generates a large charge pulse (^^lO** electrons/incident electron) which is amplified and shaped 
by a charge sensitive amplifier and one-shot circuit. The pulses are routed to the payload 
counter-board where they are accumulated and periodically read out by the payload telemetry. 

Only electron RPA's were used on this flight. Practical considerations, relating mostly tc 
the high voltage supply limitations, dictated the choice of either electron or ion RPA's but not 
both. Low energy ion measurements are effective in determining plasma flow velocities due Co 
the relatively low ratio of ion thermal speed to typical flow velocities. This amounts to an 
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independent meesuremeni of the electric held which it reitted to the drift velocity by 

PwU 

*<• ** In the ettentiully coiliiionleM plasma above 110 km. Such double measurements 

have been used to establish the credibility of double probe measurements (Carlson, 1974). 
Since the cannister included a pair of double probes, no further flow velocity measurements 
were required. It was decided that an extension of the electron measurements to lower energies 
would be preferable. 

Three RPA’s were included with aspects centered at 0*, 30*, and 180* from the payload 
spin axis. The geometric factor of each RPA was 5.7xlO"’ cm^ ster eV with a conical angular 
acceptance of half-angle 3* and an energy resolution of about 30 eV. They complete one scan 
of the electron spectrum from 6 eV to SOO eV in 128 steps each 1.024 seconds. The counts 
were collected and telemetered once every .(K)4 seconas. 


High energy elearons and ions 

Hemispherical plate electrostatic analyzers (ESA's) were employed to monitor electron 
fluxes from 350 eV to 30 keV and Ion fluxes from 470 eV to 40 keV. These detectors find 
wide application in rocket and satellite instrumentation today and no detailed description wilt be 
given here. For more information see Chase (1969). Briefly, the ESA consists of two concen- 
tric hemispherical conducting plates h< at ditTcrent potentials. Only particles of the proper 
charge sign and an energy per unit charge in the proper range are allowed to pass between the 
plates to the exit where they are detected by a CEM. 

The geometric factor of a small aperture (~ < .2) hemispherical plate ESA as shown in 
Figure II. 2b may be estimated with fair accuracy from the formula, 

S(£)- (-^'■■2/(sin<6„,„£ (11.1) 

where A is the collecting urea of the CEM placed at the exit, is the collimator acceptance 
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half angle, and E is the center acceptance energy. Note that the geometric factor increases 
linearly with the energy of the particles being analyzed. The angular acceptance in the other 

direction is assumed to be determined by the analyzer plates and is approximately 2^ at the 
extremes yielding a weighted average, for geometric factor purposes, of about The relative 

energy resolution of the instrument is essentially constant at all energies and is ~ at the 
extremes or about half that as a weighted average. 

Both ESA’s had aspects of O’with respect to the payload spin axis. The geometric factors 
of the electron and ion ESA’s were 1.8x10“^ cm^ ster (£) and 1.1x10"^ cW ster (£) respec- 
tively. Both ESA’s complete one scan over their respective energy ranges in 128 steps each 
1.024 seconds. The counts arc collected and telemetered each .004 seconds. 

Ion composition experiment 

Also included among the instrumentation was an experiment designed to yield informa- 
tion on the composition of the energetic ion fluxes in the auroral plasma. Such composition 
data can yield valuable information regarding magnetospheric acceleration and storage mechan- 
isms. For instance, multiply charged ions would be expected to have their origins in the solar 
wind with little or no storage time in the magnetosphere while singly charged ions might be 
ionospheric or magnetospheric in origin. 

The details of the instrument are to be found in Lin (1974). A similar instrument has 
been flown by Lynch et al. (1976) and detected O*, and He*^"^ ions in post-auroral breakup 
conditions. The experiment consists of a solid-state surface barrier detector sampling ions 
which have been previously sorted by a hemispherical ESA. In this case the same ESA used 
for the energetic ion measurements was used for the ion composition experiment. 

It is a well known property of solid-state detectors that the charge pulse created by the 
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pencirulion of ions with relatively low energies is a s'rong function not only of the ion energy 
but the ion mass us well. Heavy ions create relatively smaller pulses than lighter ions of the 
same energy (Potter ami Campbell. 1978; Campbell and Lin, 1973). This so culled 'pulse 
height defect" may be employed in combination with the known Elq of the ions passing 
through the ESA to yield both the charge state and muss of the ion. 

The experiment is complicated by the inherent noise properties of surface barrier detec- 
tors, especially when used to detect ions in the under 10 keV range. The intrinsic noise of such 
detectors is usually on the order of the signal resulting from penetration by a 10 keV ion. 
There are two methoils of dealing with this problem. The ions may be electrostatically pre- 
accelerated to raise their energy above the noise level or the noise level itself may be lowered 
by cooling the detector. Both of these techniques were employed in the design of this experi- 
ment. Tlie detector was lowered to a potential 5 kV below that of the exit aperture of the ESA 
thus adding an additional 5 keV’ per charge to each ion. In addition an evaporative cooling dev- 
ice was employed to cool the detector to a temperature of about -50 °C. This device was ftlled 
with liquid freon and vented shortly before cannister ejection. 

To minimire the likelihood of penetrating radiation producing spurious counts, an 
unlicorrelation detector was placed behind the main detector. The outputs of both detectors 
were routed to charge sensitive amplifiers and then to an anticoincidence circuit. 

A thrcshhold discriminator was used to eliminate response to detector noise. Those 
pulses which were ruled invalid by the anticoincidence circuit were counted and telemetered 
250 limes per second. The remaining pulses were delivered to a 12 bit (4096 bin) pulse height 
analy/cr. The I’ll A processed and telemetered 2000 events per second. Only the first event in 
the period following a readout was processed for the next readout. In addition the number of 
valid and above thrcshhold counts were telemetered 250 times per second. 
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II. B. Kleciric Held exp«rlmenls 

Three pairs of probes mounted on seif-extending booms were included on the cunnister to 
monitor DC and AC electric Helds, plasma density, and temperature. 

DC electric field 

Two pairs of the probes were designed to measure the three-dimensional electric Held. 

Referring to Figure II. I the axial booms (5 and 6) constantly monitor the component of the 

electric Held parallel to the payload spin axis while one pair of the radial booms (1 and 2) is 

» 

able to measure the two components perpendicular to this axis due to payload rotation. The 
potential difference between the pairs of probes was monitored and included in the telemetry 
stream. 

The radial booms were deployed by a centrifugally assisted feed of wire from an internal 
drum. The fully deployed length was 3.3 m tip-to-tip, but the rotation of the drum was moni- 
tored by the telemetry and less than full deployment could be accounted for. The axial probes 
deployed purely under the action of the spring loaded boom on which they were mounted, their 
fully deployed length was 2.0 m tip-to-tip. 


AC electric field 

In addition to the DC data, the signal from the electric Held probes was used to genciate 
three modes of AC Held measurements. In the first mode the signal was processed by a 
bandpass filter to roll oif frequencies less than 3 hz for the radial booms and .5 hz for the axial 
booms and frequencies greater than 1 28 hz for both pairs. By eliminating the large variations at 
the payload spin frequency, especially in the radial probes, the output signal could be greatly 
amplified. This waveform was then monitored with high time resolution (1000 samples per 
second) yielding high sensitivity to low frequency electric field fluctuations. 
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In the other two modes the signal was passed to a circuit consisting of a swept voltage 
controlled oscillator, mi.xer, and crystal Alter forming a swept narrow bandpass Alter. The log of 
the power in the passband was converted to a DC voltage and added to the data stream. Two 
such spectrum analyzers were included on the payload. The Arst swept from 0 to SO khz in SI 2 
steps with a bandpass of 300 hz completing one sweep every 2.048 seconds. The second swept 
from 0 to S Mhz in il2 steps with a bandpass of 30 khz completing one sweep every .S12 
seconds. 

ll.C. Plasma expcrlnietils y 

Two of the radial booms (3 and 4) similar to those used for the electric Aeld experiments 
were employed as Langmuir probes to monitor plasma temperature and density. For a detailed 
discussion of Langmuir probe operation see Lochte-Hoitgreven (1968). BrieAy the current- 
voltage characteristics of the probe when immersed in a plasma are determined primarily by the 
electron density and energy distribution. By controlling and varying either the voltage of the 
probe or the current it draws from the plasma, and monitoring the other an experimental 
characteristic curve may be obtained which yields both the temperature and the density. 

The proper operation of a spacecraft mounted Langmuir probe depends on the assumption 
that the payload itself is well "grounded" with respect to the plasma. In practice this means that 
the conducting surface area of the payload must be much larger than that of the probe to insure 
a low impedance return path to the plasma for the current collected by the probe. Failing this 
the spacecraft itself, which ordinarily develops a sheath potential on the order of the electron 
thermal energy, may be driven to much larger voltages. For example, if the probe is held at a 
large positive voltage it collects essentially all of the electrons and none of the protons in its 
vicinity. This current must be balanced by an equal current of ions to or photoelectrons from 
the payload surface. Due to the much slower thermal speeds of ions, cither a much larger sur- 
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face area or a very large negative voltage is required of the payload. In fact in the very attempt 
to drive the probe highly positive, the payload may be driven highly negative while the probe 
remains near the plasma potential. Spacecraft charging problems cause difficulties with low 
energy particle detectors which depend on the spacecraft to maintain its potential close to that 
of the ambient plasma. As will be seen in the next section this appears to have happened dur* 
mg the night of the probes described here with the result that the RPA's did not function prop- 
erly at energies above 100 eV. 


By applying a positive voltage well in excess of the anticipated electron thermal energy to 
a Langmuir probe, the current /q, referred to as the electron saturation current, is collected. 
For a spherical probe in the ionosphere, where the Debye length is a fraction of a centimeter 
and sheaths arc thus very thin, /q is given by 


/o - 


2KT, 


ftu 


( 11 . 2 ) 


where A is the piobe area and /w,. is the electron mass. In this case the probe was held 6.2 
volts positive with respect to the payload. This is well above the thermal energy of the ionos- 
pheric electrons which is usually a few tenths of an eV. The current collected by this probe is 
converted to a voltage and processed by a log ampliher. Both a direct output sampled 1000 
times per second and an AC coupled high gain output sampled 2000 times per second are 
included in the telemetered data. 


The current collected by this probe is used to control a constant current source which 
forces the other probe to collect a small fraction of the electron saturation current, /’. This 
probe operates in the exponential region of the Langmuir probe characteristic curve. In this 
region 
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A small sinusoidal variation of amplitude .1/' is added to the bias current at a frequency 
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- I kHz. This produces a variation in the voltage of the probe given by 

KT 

5J/-.I — ^sino)/ (II.4) 

e 

This voltage variation is amplihed, demodulated, and included in the telemetry data with a sam- 
ple rate of 2S0 per second along with an AC coupled high-gain channel of the same data sam- 
pled SOO times per second. Using Equation II. 2 and Equation II.4 these data determines the 
electron temperature and density. 


II.D. Magnetic field experiment 

A two axis flux-gate magnetometer was included on the cannister primarily for attitude 
analysis. The axes coincided with the x and z axes as shown in Figure II. 1 The DC voltage 
output of each axis was converted by a 12 bit analog to digital converter and telemetered 250 
times per second. The instrument had a sensitivity of 10 nT but was limited by the range of 
the analog to digital converter to a resolution of one bit which corresponded to 28 nT. As such 
it was potentially useful in detecting the magnetic perturbations due to ionospheric sheet 
currents which often are several hundred to a thousand nT. 



III. LAUNCHES 


In February and March of 1976 a campaign was mounted to launch the hrst two Porcu- 
pine payloads from Kiruna, Sweden. Porcupine I was launched during this campaign. Unfor- 
tunately a malfunction of the Aries launch vehicle resulted in the destruction of the payload 
about one minute after launch and*no useful scientific data were returned. The planned second 
launch was postponed pending study of the Aries vehicle. 

The following year, tests revealed the probable cause of the malfunction - inadequate heat 
sinking on the motor dome. The problem was corrected and a second campaign was mounted 
in February and March of 1977 to launch the second payload. 

On March 20, 1977 at 19:22 UT (about 22:30 magnetic local time) Porcupine 2 was 
launched into a quiet arc. The arc diminished in strength in the minutes prior to launch and 
was too faint to be reliably recorded on all-sky cameras during the flight. Thus no orientation 
data are available from this source. Poleward of the rocket position, however, a visual form 
was observed which ran from somewhat south of west to somewhat north of east. 

Both UCB instrumented cannisters were ejected properly but for a still undetermined rea- 
son data were received from only one of the pair. In addition, boom deployment problems 
resulted in the re-orientation of this cunnister’s attitude from its planned configuration. 

The cunnister was programmed to deploy the four radial booms first to enhance the rota- 
tional stability about the intended spin axis with the two axial booms deploying some 20 
seconds later. Apparently the r.idial boom deployment mechanism locked up preventing full 
extension and possibly causing unsymmetrical extension. When the axial booms were released 
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the cannisier rapidly entered a period of unstable gyration. After about 40 seconds it settled 
into a very stable flat spin about an axis perpendicular to the axis of symmetry (the 2 -axis as 
shown in Figure II. 1) and remained in this attitude for the rest of the flight. A good attitude 
analysis is thus available for times later than TAL (time after launch) 180 seconds. No 

attitude analysis is available prior to this time. 

A second problem was the apparent charging of the payload. As was noted in the previ- 
ous section, the potential of the spacecraft is determined by the requirements of current balance 
between its body and the plasma. The usual method of insuring that the spacecraft body 
remains at ambient plasma potential is to make the area of the spacecraft skin much larger than 
the ratio of the electron to ion thermal speed times the probe surface area. In the case of the 
ejectable probes this condition was no more than marginally met. If the current can not be bal- 
anced with the spacecraft at plasma potential it will be driven to lower potentials until the 
reduced electron current to the probe and the enhanced ion current to the spacecraft body 
satisfy the current balance condition. 

As a result of this charging, large fluxes of thermal ions were allowed past the small elec- 
trostatic barrier which would otherwise screen thermal ions from entry to the retarding potential 
analyzers (see Figure II. 2a). Once inside the focusing region of the detector, they were drawn 
to the highly negative walls emitting secondary electrons. These were then drawn to the -f ISOO 
Volt entrance of the CEM producing spurious counts. These effects were most noticeable dur- 
ing the high energy end of the RPA sweep because it was at these times that the walls of the 
f(K'Using region were at the largest negative potentials and thus, that the collisional energies of 
the ions were the greatest. The result is that low energy electron data were lost at energies 
above 100 eV. 

The tr;gectory of the cannister is determined from radar data obtained for the main pay- 
load. Since the physical separation of the cannisters from the main payload was never more 
than two km this assumption causes no problems. The payload reached a maximum altitude of 
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459 km at f ■■ 370 seconds. At r ■■ 373 seconds the barium shaped charge experiment was 
detonated. Reentry and loss of signal occurred at about t - 660 seconds. 

The tnuectory is shown in Figure lll.l. Note the exaggeration of the horizontal scale. 
Also included on this flgure are lines of constant L or magnetic latitude. Due to the extremely 
small horizontal range of the flight, the post -apogee portion of the flight provided a magnetic 
retrace of the pre-apogee portion. This fact was useful in determining the time scale for 
changes in the observed auroral precipitation. 

In Figure III. 2 the 110 km altitude magnetic projection of the trajectory is displayed. The 
dashed line represents the orientation and equatorward edge of the arc encountered during the 
flight. The orientation is determined from theoretical considerations os explained in the 
analysis section and corresponds to the approximate orientation of the visual form observed 
poleward of the rocket position. 
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IV. DATA 


IV. A. Overview 

Electron measurements 

Figure IV. I serves us a broad introduction to the important features in the flight data. In 
all three panels the horizontal coordinate is flight time and runs from r » 120 seconds (shortly 
after cannister ejection) to t - 660 seconds (re-entry), The top two panels display the electron 
dilTerentiul energy flux. 4>, as a function of electron energy and time. The top panel is com- 
posed of data from the electron ESA while the middle panel is from one of the RPA's. The 
RPA data are cut off at 100 eV due to the malfunction referred to in the previous section. The 
bottom panel displays data from the high frequency electric field spectrum analyzer. Darker 
regions correspond to greater intensities in all three panels. 

The conventional parameter used to report directional particle fluxes is J, the dilferential 
particle flux, measured in particles /cm* s ster eV. Since <1> -• EjiE), the dilTerential energy 
flux tends to make the secondary peak in the electron spectrum stand out more clearly. The 
auroral fluxes arc clearly defined in the ESA data with the spectral peak reaching a maximum of 
about 4 keV at t - 420 seconds. In fact, except for a brief period from 160 to 230 seconds, the 
payload was constantly within u region of energetic electron precipitation. The energetic fluxes 
display a periodic fluctuation in intensity which is due to a beating effect between the energy 
sweep frequency of the detector itself and the payload rotation frequency. The detector had a 
view angle of 90° with respect to the final spin axis and the spin axis formed an angle of 30° 
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with the magnetic Held. Thua during the course of one rotation the detector umpled the parti* 
cle flux at pitch angles from 60* to 120*. As wilt be shown later in the data presentation, the 
electron fluxes are greater over the downward pitch angles. This results in apparent temporal 
fluctuations due to the variations in detector view angle during successive sweeps. The pilch 
angle distribution of the energetic electrons will be discussed later in this section. 

The RPA data essentially follow the ESA data displaying higher intensities in the low 
energy electrons corresponding to increases in the primary beam energy. This result is 
expected due to the enhanced production of secondary electrons through collisions with neutral 
particles in the ionosphere. These secondary electrons are much more isotropic and thus no 
beating effect as seen in the ESA data occurs. 

Between .160 and 420 seconds there are four brief intensifications noticeable in the RPA 
data. The first of these is an artifact of the barium shaped charge experiment which was 
detonated at t - 373 seconds. The apparent intensification at about 40 eV is merely a result of 
the fact that the shock front of the expanding barium cloud happened to pass the location of 
the cannister at a time when the RPA was sampling fluxes at that energy. The other three 
events, all at an energy of about 2S eV, are real. In the first and last of these events the 
intensification lasts for no more than the duration of a single energy sweep (% one second) and 
thus we cannot be sure that we are not observing a very brief enhancement that is occurring at 
all energies. However in the second event the intensification exists to some extent for at least 
two sweeps. The two other RPA*s not shown here show similar enhancements at the same 
times and over a variety of pitch angles. This evidence and the fact that all the events occur at 
similar energies seem to imply that the phenomena are similar and do represent a brief 
enhancement in the electron flux over a limited energy range near 25 eV. 
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AC Electric fieM 

The bottom panel of Figure IV. I displays the power in the high frequency electric Held 
Huc;uations as observed by the spectrum analyzer. There are several features of this plot which 
require explanation. The periodic vertical sets of five dark rectangles are frequency calibration 
markers generated by an omniard oscillator to provide reference frequencies for inflight calibra- 
tion. The dark bands which run across the top of the panel are short wave radio transmissions. 
Their momentary disappearance during re-entry is due to the elevated cutoff frequency pro- 
duced by the enhanced electron density in the tower ionosphere. Also plotted on this panel are 
the electron plasma frequency as determined from electron density measurements and twice the 
electron gyrofrcqucncy as determined from the model magvunic field, both at the payload posi- 
tion. These frequencies will be important in the analysis of the spectrum analyzer data. The 
cyclotron frequency is observed to change slowly as the payload moves in altitude. Likewise 
the electron plasma frequency decreases as the payload moves to regions of lower electron den- 
sity at higher altitudes. The major perturbation in the middle of the flight is once again a result 
of the barium shaped charge experiment. The initial drop in plasma density is a result of elec- 
tron attachment to the chemical byproducts of the explosion. Within a few seconds photoioni- 
zation of the neutral barium produces a dramatic enhancement in the plasma density which 
then decays slowly over the course of the next 40 seconds. 

There are several instances of isolated bursts at the plasma frequency and twice the elec- 
tron gyrofrcqucncy in the early part of the flight. By far the most dramatic event occurs after 
the barium detonation around t - 400 seconds. Large amplitude waves occurred along the 
line for about IS seconds following the decay of the locally elevated plasma frequency to a 
value below twice the electron cyclotron frequency. These emissions culminated in a brief 
burst of waves over a wide range of the spectrum. We will return to these data and examine 
them in much greater detail later. 
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DC Etednc FieM 

Figure IV. 2 displays the southward and westward components of the electric field perpen* 
dicular to the magnetic field and the total electric field as a function of time as observed in a 
geographic (corotating) frame of reference. The observed component of the electric field paral- 
lel to the magnetic field is consistent with zero. No data are presented prior to r ■■ ISO seconds 
due to the lack of an adequate attitude solution during the early portion of the flight. The large 
perturbation in the data at t - 373 seconds is once again an effect of the barium shaped charge 
experiment which produces a momentary perturbation of the plasma flow velocity and thus the 
observed electric field. 

Figure IV. 3 shows the magnetic projection of the observed DC electric field to an altitude 
of 110 km. The plasma flow implied by the resulting ExB drift is generally equatorward and of 
much larger magnitude outside the arc than inside. 

/on measuremems 

Throughout the flight the counting rates detected by the ion ESA were too far too low to 
be statistically significant without averaging over long times. Even then most counts arc due to 
detection of scattered ultraviolet from those times that the instrument viewed angles close to 
the sun. Figure IV. 4 is an example of a spectrum obtained by averaging ion data for 50 seconds 
starting at t “ 250 seconds. This spectrum is representative of that seen throughout the flight. 
As a result of the general luck of energetic ions, the ion composition evperirnent also yielded 
no data of interest. 


Maitnvtomeier data 

As was discussed in the instrumentation section, the magnetometer data were used pri- 
marily for attitude analysis. However, by subtracting the model magnetic field from the 
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obierved data, ftcld-uligned and ionospheric sheet currents may be revealed as perturbations in 
the residual held. The residual noise limited the resolution of magnetic perturbations to about 
25 y. This would correspond to a sheet current of about .02 A/m. No fleld*aligned sheet 
currents of this magnitude were discerned in the data. Just before the end of the data, a per- 
turbation of about 90 y was observed. This is attributed to an E-layer sheet current of about 
.08 A/m which wus traversed on reentry. 


IV. B. Antlcorrelation of electric Held and electron precipitation 

Comparison of the top panel of Figure IV. 1 and the bottom panel of Figure IV. 2 snows 
that the variations in the total electric Held closely mimic, in an inverse sense, the variations in 
the peak energy of the auroral electrons. In Figure IV. 5 the inverse of the DC electric field is 
plotted along with the secondary peak energy of the simultaneously observed precipitating elec- 
trons. It has been noted in numerous sets of auroral data that the ionospheric electric field 
dccrca.ses and sometimes increases (Dc La Beaujardierc, 1977) within auroral arcs but it is 
interesting to note the remarkable degree of anticorrelation which occurs within this arc. The 
two temporary retreats in the peak electron energy as the aurora gains in strength are repro- 
duced quite faithfully as temporary recoveries in the electric field strength which generally 
decrciises during this portion of the fiight. The detailed relationship between the particle fluxes 
and the electric field will be explored in detail in the analysis section. 


IV.C. Electron distribution function 

One of the more interesting results of this flight wus the observation of changes in the 
distribution function of the precipitating electrons apparently associated with the occurrence of 
the barium plasma injection. In the statistical description of plasma dynamics the distribution 
function, /(r.v). which is the number of particles per unit volume in phase space, is used. 

r 
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Typically, panicle flux measu 1*60161111 are reported in terms of the differential particle flux, 
J(E,a), which is the number of particles hitting a unit area per unit incident solid angle per 
unit time per unit energy at the energy E and the pitch angle a. With the appropriate change 
in variable dependences it is easily shown that these two functions are reiated by 



In the discussion which follows we will use the quantity JIE when referring to the distribution 
function. 

Due to the unanticipated change in attitude of the cannister, the electron ESA monitored 
electron fluxes at pitch angles between 60 and 120 degrees. Figure IV. 6 displays two sets of 
electron distribution functions as measured by the electron ESA, one from before the barium 
injection and one from afterwards. These data are obtained as the cannister rotates and no 
attempt to sort them by pitch angle has been made here. Nevertheless a clear difference is 
noticeable in the fact that the pre-barium distribution functions display a marked peak at 2 to 3 
keV while th' post-barium distribution functions do not. The one non-peaking case in tho, 
pre-barium examples occurs at a time when the payload attitude gave the electron ESA a view 
of electrons with pitch angles greater than 90 degrees, that is electrons moving upward. 

These data may be more clearly presented in the form of a polar plot of iso-intensity dis- 
tribution function contcjrs as in Figure IV. 7. This figure displays distribution function con- 
tours in electron velocity space. The left (right) half of the figure represents the distribution 
function before (after) the barium injection. Notice that data are obtained only between 60 and 
1 20 degrees. The rest of the figure is drawn in as an educated guess based on an extrapolation 
of these data and results of previous observations and should not be taken too seriously. On 
this plot the thermal electrons are represented by an intense narrow spike at the origin while 
the secondary electrons formed by the collisional interaction of the primary beam with the 
atmosphere is seen as the broad, downward sloping portion at moderate velocities. 



It lltM 


27 






The prc-burium distribution function contours ure characterized by a ridge ut downward 
pitch angles which is eroded into a plateau following the barium injection. The ridge begins to 
recover about 60 seconds after the barium detonation but never regains its pre-barium magni- 
tude. There is good evidence that the distribution actually peaks close to 90 degrees. 

Such secondary peaks in particle distribution functions are well known sources of free 
energy for the growth of plasma waves. The deterioration of the electron distribution peak wilt 
be explored in more detail in the analysis section. 

IV. D. Iligh-frequency plasma waves 

Figure IV. 8 displays twelve succc.ssive half second sweeps of the high frequency spectrum 
analyzer including the period of intense emissions near t - 400 seconds referred to earlier. In 
the first sweep the most prominent features are a sharp line at the upper hybrid frequency (= 
2.4 Mllz) and . broader line at the plasma frequency (= 2.05 Mllz). In successive sweeps the 
upper hybrid emission grows to very large magnitude and then dies away. Figure IV. 9 is u 
detailed version of the se. enth of these sweeps and shows the richness of the observed spec- 
trum very clearly. 

The upper hybrid line is so sharp that it is not resolved by the instrument. Assuming its 
width to be equal to this resolution, a lower bound for the line strength may be computed. 
Such a calibration is used for the abscissa in Figure IV. 8 and Figure IV. 9 instead of the usual 
mV/m /Hz'-. A line strength reading of I mV/m corresponds to a spectral density of .033 
mV/m /kHz '. The upper hybrid emission line strength calculated in this manner is greater 
than 20 mV/m during the strongest portion of the emissions. 

As the upper hybrid emission grows waves at a large variety of frequencies are observed. 
At least two and perhaps three sidebands emerge at 200 kHz intervals on cither side. On the 
low frequency side these sidebands ride on the back of a broad plasma frequency emission. A 
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pair or lower frequency emissions at 200 and 400 kHz. several small very high frequency (4-S 
MHz) oscillations, and Anally a large emission feature at = 650 kHz are also observed during 
this period. 

Several of these emissions At well into the picture of a parametric decay process taking 
energy from the upper hybrid pump wave and distributing it into the sidebands and the 200 and 
400 kHz emissions. Also the emissions appear to be stimulated in some manner by the pres- 
ence of the barium cloud. These waves, their relation to the barium cloud, and the changes 
observed in the electron distribution function will be analyzed more thoroughly in the analysis 
section. 

IV. E. Summary 

To review the main points of the data, the payload was in a soft precipitation region equa- 
torward of a quiet arc from r - 165 seconds to t - 230 seconds. For the rest of the flight it 
encountered auroral electron precipitation at energies up to 4 keV. The DC electric fleld mag- 
nitude showed a remarkable degree of anticorrclation with fluctuations in the peak energy of 
the precipitating electrons. At t - 373 seconds a barium shaped charge experiment was 
detonated causing a large plasma density enhancement surrounding the payload. Immediately 
after this the local electron distribution function lost the prominent peak which it had previ- 
ously displayed at downward pitch angles. Several brief enhancements in the flux of 25 eV 
electrons were observed. A strong emission at the upper hybrid frequency began when the 
local plasma frequency returned to the low side of the second electron cyclotron harmonic and 
ended in an intense burst of plasma emissions at a large number of frequencies. 


In the next section we will discuss possible mechanisms for these observations. 



V. ANALYSIS 


V.A. Oil the control of the ionospheric electric Held 

A simple model and some considerations 

The remurkuble anticorrelation between particle precipitation and the electric field seems 
to imply a controlling influence of one phenomenon over the other. The relationship at once 
suggests the possibility that the ionospheric conductivities are modulated in a temporally 
responsive fashion by the instantaneous electron precipitation and its ability to enhance the 
plasma density in the lower ionosphere. The presence of a large scale ionospheric current sys- 
tem might then determine the ionospheric electric held and produce the type of inverse rela- 
tionship observed here. 

What might be required to produce the observed aniicorrelation? One could at first be 
led to consider the relationship of the ionospheric electric fleld to the field-aligned electiic fleld 
which proouces the energetic electrons. An application of J*E r/l - 0 to closed paths linking 

t 

the magnetosphere to the ionosphere might be expected to yield a type of inverse relationship 
between the field-aligned potential drop and the ionospheric field. As a result, a relationship of 
the form E + constant' (PeakEnergy) - constant might be obtained. A quantitative model of 
this sort is impossible due to our lack of knowledge of the boundary conditions, namely - the 
magnetospheric potential configuration. However, for the simple assumption that E,„„^ - 0 or 
- constant. Figure V.la and b displays the resulting ionospheric and peak electron energy 
profiles. Figure V.lb differs from Figure V.la only by virtue of an offset imposed on the 
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ionospheric electric held by the magnetospheric potential gradient. No inverse relationship is 
evident in either of these cases. On the other hand, Figure V.lc shows that it is possible to 
construct a potential configuration for the magnetosphere which produces the postulated rela- 
tionship (or any other relationship). We view this possibility as remote, especially when the 
contortions required to reproduce the complex peak energy profile of this arc are considered. 

An alternative model could be constructed as follows. First the ionospheric conductivities 
should vary in direct proportion with the peak energy of the precipitating eiectrons. This would 
be due to the enhanced ionospheric electron density created by the energetic electron precipita- 
tion. Second the ionospheric currents perpendicular to B should be constant in both time and 
space. That is, the currents are determined by large scale features of the magnetosphere- 
ionosphere interaction and are relatively unaffected by the arc. Implicit in this assumption is 
the requirement that negligible field-aligned currents flow in the region of interest since they 
would lead to a divergence of and therefore a variation in the ionospheric current system. Even 
though substantial Acid-aligned currents are known to flow in auroral arcs, this assumption is 
not unreasonable. Maximum Aeld-aligned current densities are typically a few juA/m^ and even 
a 10 km wide region of such a current density would produce only a few hundredths of an A/m 
of field-aligned sheet current. Ionospheric sheet currents can easily reach an A/m and thus will 
be relatively unaffected by such a Aeld-aligned sheet current. Under these assumptions the 
electric Acid would be inversely proportional to the conductivity and, thus, the peak energy of 
the precipitating particles. 

Electric Acid measurements in the ionosphere are highly frame dependent and thus any 
attempt to link measured electric Acids with other measurements must be approached with cau- 
tion. For calculations involving the classical Hall and Pedersen conductivities the proper frame 
in which to work is that of the neutral particles. Unfortunately this frame is neither well- 
dcAned, due to large shears in the neutral wind velocity with altitude, nor even measured in 
most cases. Thus the effective electric Aeld for current calculation purposes may be altitude 
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dependent due only to the altitude dependence of the neutral frame itself and, worse yet, com- 
pletely unknown. 

Furthermore there are several problems with this model. Time variations in the precipita- 
tion arc not considered. Also the fact that the current is actually related to the electric field 
through a tensor and not a scalar conductivity compli'^ates the problem. 

When the electron precipitation varies there is a time lag associated with the build up or 
decay of the plasma density. The equation governing the electron density, », in the presence 

of a production rate, </, and a recombination coclUcient, a , is 

»• 

~ m, q — (V.l) 

dl 

The recombination rate is elTectivcIy a function of altitude through its dependence on the tem- 
perature. The production rate is also a function of altitude. Due to the non-linearity of Oqua- 
tion V.l the time scales for build-up and decay of the electron density are functions of the 
strength of the aurora as well as the altitude. For stronger precipitation the recombination term 
becomes important more quickly and the electron density "tracks" the changes in the precipita- 
tion more closely. Jones and Rees (1973) used a model ionosphere and a typical auroral pre- 
cipitation profile to calculate the time constants for electron density build-up and decay. The 
times vary from 5 to 70 seconds in the altitude range considered. The fastest time constants 
occur in the E-region which is where the miyor contribution to the ionospheric height 
integrated conductivities is found. Thus one might expect typical lags on the order of ten 
seconds between precipitation changes and conductivity changes. Note that this effect is impor- 
tant even in the case of steady precipitation if there is motion of the ionized component relative 
to the arc. For example, How across an arc would be expected to produce enhanced conductivi- 
ties downstream from the arc and anomalously low conductivities within the arc. Since plasma 
flow speeds may be as high a« ! or 2 km/s and auroral thicknesses can be on the order of a few 
km, the plasma may be able to completely traverse an arc in a time less than the build-up or 


V. 


32 






decay time constant. This would lead to reduced conductivities within the arc and an extended 
downstream region of enhanced conductivities relative to those expected in a non*row situa- 
tion. 

Another complicating effect of time variations is seen in the fact that the electric field 
itself may not be able to change instantaneously even if the conductivities could. In the 
absence of wave turbulence generated anomalous resistivity along magnetic field lines, the 
ionospheric electric field is coupled to the magnetosphere due to the high fleld-aligned D.C. 
conductivity of the magnetospheric plasma. Changing the ionospheric electric field is then 
equivalent to altering the plasma convection velocity and thus the momentum of all the magne- 
tospheric plasma contained in the flux tube linking the region of the ionosphere under con- 
sideration. Anomalous resistivity effectively decouples the motion of the magnetospheric 
plasma from that of the ionospheric plasma by enabling field-aligned potential drops to exist. In 
this case it might appear that rapid changes in the ionospheric electric field would be possible 
since the inertia of the magnetospheric plasma has been removed. We must remember, how- 
ever, that we are talking about time variations in the ionospheric electric field and thus we need 
to consider the A.C. field-aligned conductivity. This A.C. conductivity is determined by the 
properties of Alfven waves which transmit the perturbation electric field along the field lines. 
Thus, for time varying electric fields, the inertia of the magnetospheric plasma still plays a role. 
When the precipitation and thus the ionospheric conductivities change in a limited spatial 
region, a non-zero divergence of the horizontal current system in the ionosphere creates charge 
layers at the boundaries. This divergence cannot persist and is eventually eliminated in one or 
a combination of two ways. If A.C. field-aligned currents are somehow prevented from flowing 
the ionosphere will polarize until the resulting electric field once again permits a divergenceless 
horizontal current system. If, on the other hand, A.C. field-aligned currents are completely 
unimpeded they might account for any divergence with no need for polarization. 
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In fact neither of these mechanisms can exist in the absence of the other. In the Arst case 
complete polarization would imply a decoupling of the ionospheric electric held from that of the 
magnetosphere. This situation would soon be eliminated by Aeid-aligned currents driven by the 
resulting Aeid-aligned potential diAference. The polarization Aeld created in the ionosphere will 
eventually be communicated to the magnetosphere. In the second case, Aeld-aligncd currents 
imply the Aow of charge which will inevitably be accompanied by a perturbation electric Aeld. 
In other words the ionosphere must polarize to some extent to drive the Acid-aligned currents. 

Thus the only self-consistent solution of the problem is a combination of the two eflccts. 
The question of which is the dominant process can be answered by looking at the relative 
impedances of the two return paths available for the polarization charge created by the initial 
divergence. Mallinckrodi and Carlson (1978) have shown that the conductivity for Aeld-aligned 
currents associated with perturbation Adds, the so-called wave conductivity, is on the order of 
.5 mho while auroral ionospheric conductivities are typically 10 mho. Thus, in fact, A.C. Aeld- 
aligncd currents are also relatively highly impeded. Over the course of time a new D.C. 
conAguration is established. Goertz and Boswell (1979) have determined that the steady state 
currents arc iichicved on time scales of the order of hours. As a result the ionosphere docs 
polarize but the inductive effects of the momentum of the plasma on the magnctosphcric Aux 
tube delay the establishment of the steady state electric field. Once again this time dependence 
need only occur in the frame of the plasma and so convection plays an important role in this 
process as well. 

Finally the model is incomplete because of its failure to account for Hall currents which 
can be even greater than the normal Pedersen currents. Since the ionospheric currents do not 
flow in the direction of E and since variations in the auroral electron density affect the Hall and 
Pedersen conductivities differently, one cannot simply assume that a constant and uniform 
ionospheric current traversing ionospheric plasma density variations will lead to an inversely 
related electric field. Obviously we must be more precise in the mathematical description since 
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the relationship! are more complicated than indicated by our simple model. 

The numerical simulation and application to flight data 

In an attempt to understand the observed anticorrelation, a numerical procedure was 
developed to simulate the electrodynamics of auroral ionospheric phenomena. The routines are 
described in detail in the appendix and the reader is referred there for a more detailed discus- 
sion of the calculational procedure. In the following discussion we refer to the coordinate sys- 
tem defined there and shown in Figure V.2. The assumptions of the numerical routine limit us 
to considering only those cases in ’which there is a sensible "arc frame" in which all quantities 

A 

are time independent and in which ~ (the longitudinal gradient) vanishes. Both of these 

ay 

assumptions are fairly well justified due to the fact that the visible aurora displayed the typical 
longitudinal extension and lack of dynamics of a quiet arc. 

Under these assumptions and the further assumption of a negligible held-aligncd com- 
ponent of the electric field (well justified below 2S0 km), it is required that be constant 
throughout the region and time of interest. The electric field data from the flight indicate a 
fairly constant component in a direction approximately 30° south of west. Figure V. 3 displays 
the electric field of Figure IV. 2 rotated by 30° into this orientation. Notice that, with the 
exception of the barium perturbation, £,. is fairly constant for a good portion of the flight. The 
particular arc that the rocket penetrated was too faint to be observed by all sky cameras. How- 
ever this direction corresponds well with the observed orientation of the arc poleward of the 
rocket position. Since multiple arcs generally align themselves in similar directions we choose 
this orientation as the most likely direction of extension for the arc. 

Next we must determine the motion of the arc frame relative to the geographic frame in 
which the electric field measurements are reported. Figure IV. 1 reveals a certain degree of 
symmetry in the particle precipitaiion about / - 420 seconds. This time is very close to the 
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time at which the payload reverses its direction magnetically as shown in Figure III. I and Figure 
III.2. Thus a likely possibility is that the arc is essentially motionless in a geographic frame and 
the rocket merely penetrates to a position which may even be equatorward of the maximum 
energy electron precipitation. Notice that only motions in the x-direction are perceptible due to 
the assumption of vanishing y dependence. Thus the y-component of the arc’s motion is arbi- 
trary at this point. Some time variation clearly occurs since the particle spectra are slightly 
different during the downward portion of the flight when the payload retraces the magnetic 
coordinates of the upward leg. These differences are slight and occur on a time scale of 
minutes. As was pointed out above, the electron density in the ionosphere requires times of 
the order of tens of seconds to adjust to changes in ionization rates. Thus the steady state 
assumption seems well justified. 

The data between 220 and 410 seconds were used to obtain average precipitating particle 
spectra at one kilometer intervals from the equatorward edge of the arc to the maximum pene- 
tration position. This corresponds to a total distance of 32 kilometers perpendicular to B. 
Notice that during the interval between 380 and 44S seconds the payload's motion transverse to 
B was limited to about one kilometer. 

Using these spectra and an ionospheric model, ionospheric production rates, </, were cal- 
culated. The program uses the collisional ionization algorithm of Rees (1963) in a manner 
similar to Evans (1977). Although there is some evidence that r.on-collisional processes may 
be important in limiting the altitudinal extent of the ionization produced (Stenbaek-Nielsen and 
ilallinan, 1979), such processes seem to be more important in pulsating auroras. Once again, 
see ,nc appendix for more details on the calculational procedure. To calculate the equilibrium 
electron densities, we must know v,, the velocity of the ionized component relative to the arc 
frame. The y-component of v, is irrelevant for the calculation of the convective derivative of 
the electron density so we choose it to be zero and let r, be the component of the ionization 
drift in the -I-a- - direction. As a first approximation we will assume no drift. In this case there 


is no problem with time constants for ionization build*up and decay and the steady state may be 
assumed. In the steady sute, the electron densities are determined from 

n >/fl7o (V.2) 

These densities are then used to calculate the height integrated Hall and Pedersen conductivi* 

ties. Figure V.4 displays the calculated height integrated conductivities as functions of position 

within the arc. 

Now note that in our model 

A/. - -J. (V.3) 

where 1 is the height integrated current perpendicular to B and is the field aligned current 

A 

density into the ionosphere. (Recall that - 0 in this model.) In accordance with the mag' 

0>' 

netometer data which indicate no appreciable field-aligned currents (see discussion of magne* 
tomeier data in the daia section), we set = 0 yielding A constant. Notice that A need not 
be constant with respect to x. 

The use of our conductivities to determine the relationship between the electric field and 
the particles is hindered by one major gap in our information - the neutral wind speed. The 
relationship 



is true only in the rest frame of the neutral particles. To perform the calculation in the arc 
frame we must correct the electric field and obtain 

I-£(E+— xB) (V.5) 

c 

where v„ is the neutral wind velocity with respect to the arc frame. Note that the current is 
independent of reference frame under the assumption of quasi-neutrality for the ionospheric 


plasma. 
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Notice ihui in Equation V.5 v« is usumed to be independent of altitude. This is clearly 
not generally the case, however the fact that the ionospheric current flows predominantly in an 
aliiiudinally thin layer of about 20 km around ihe position where the conductivities peak might 
mean that as tong as an average neutral wind velocity in that region is used, the error could be 
small. In fact even in this limited region the neutral wind shears can be as large as 100 m/s 
(Roble el al., 1979) but with no data available we must make the assumption as noted. We 
further assume v„ to be independent of.v. Then from Equation V.3 and Equation V.5 we And 
that specifying at one position determines f, al all other positions. Speciflcally, 


£-,'(0)I,.(0) - I/,(0) - Z„(x) 

E:(x) (V.6) 

Zp(x) 

y 

where E'-E+~^xB It is convenient al this point to specify the ^'-component of the the neutral 

wind velocity relative to the arc frame to be equal to zero. This can be done without loss of 
generality. Figure V.S shows the relative velocities between Ihe geographic, arc, and neutral 
frames. Since the electric held measurements are presented in a geographic frame, we must 
transform them to the arc frame. With this choice of neutral wind speed, which is equivalent 
to specifying the as yet unspecified .v-componcni of the arc frame velocity relative to the geo* 
graphic frame, it is a simple matter to transform the geographic electric field to the arc frame. 
The transformation is given by 


£, - £,'' + —B 

c 

£, - E! ' - —B 
c 


(V.7) 


(V.8) 


where </, is the .t'*componcnt of the neutral wind velocity relative to the geographic frame and 
I/, is the .V -component of the arc velocity relative to the geographic frame. 

With this choice Equation V.6 can be written 


£,(.v) - £,(0)|^-(£, - --£) 


I/,(0) - 


Zp(x) 


(V.9) 


f 


38 

where it lo be uken positive in the <f jr>tiirection and neiative in the •jr>direciion. In our 
Art! ipproximstion we have already taken v. * ■■ 0. We will now alto attume - iv “ 0. 

Thus we ctn use the electric Aeld data of Figure V.3 . the conductivities of Figure V.4 , and 
Equation V.9 to compare the model to the data. To do this the position of maximum poleward 
penetration at r ■■ 410 seconds is ciu ien as jt -• 0. The .v ‘component of the electric Aeld at 
this time was esseniialiy >:ero so we set (0) ~ 0. Also, in the time interval of interest E^ was 
relatively constant at abou^ 19 mV/ni »o we choo«e this value for The results of applying 
Equation V,9 to these values are shown in Figure V.6. in this Agure £, (a) as calculated is 
comp>red with the observed E,(x). The transformation of the data of Figure V.3 to that of 
this Agure involves the use of only that portion of the data between i >• 220 and 410 seconds, 
reversing its order, and accounting for the fact that the jr>component of the payload velocity 
was not constant during this interval. 

As expected, the Agure shows general agreement in the broad features but there is room 
for improvement. In particular the Acid equatorward of the arc is much too great. Also the 
data curve is much smoother than ‘ calculated curve. 

Both of these problems can be overcome by recognizing the effects of Anite arc motion 
and plasma convection. As noted previously, these have the effect of introducing time varia- 
tion into the precipitation seen by a given element of the ionospheric plasma. Then build-up 
and decay times become important and the plasma density at any point is a function of the pre- 
cipitation it has experienced over the past minute or so. Thus we now need to consider more 
realistically the probable values for v,, v„, u^, and w,.. We will start by pointing out some of the 
general considerations which limit the possible values of these parameters. 

First, neutral wind speeds rarely exceed 200 m/s (Hays ei ai. 1979). As a result m, and 
t (/, -t- v„, the components of the neutral wind velocity in a geographic frame, should be no 

larger than this. Notice that itself is not restricted since it is the speed of the neutral winds 


relative to the arc frame. 
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Second m« is probably not very large. As was pointed out before, the electron precipita- 
tion data imply that the arc was both fairly steady and motionless in the geographic frame. 
Assuming as we have that the variations in the data are spatial and not temporal, an arc *’pecd 
of more than about 100 m/s would have shifted the observed symmetry noticeably iiway from 
the point of maximum poleward magnetic position. 

Finally, what can we say about the at - component of the ionization drift velocity in the 
arc frame? At high altitudes v, is clearly just the normal I'xB plasma drift. At very low alti- 
tudes it is the neutral wind speed since the plasma becomes collision dominated. What about in 
between? At v/hat height does the transition take place and how does •>, behave in the transi- 
tion region? As discussed in the appendix, the transition takes place belt w 100 km. Since 
almost no ionization occurs below this altitude in an arc as weak as this one it is reasonable to 
assume v, to be equal to the .v-component of the plasma drift in the arc frame. 

For the moment let us consider the effect of this non-zero v, on the previous results. 
Since - 19 mV/m. we actually have a southward r, of about 380 m/s. At this speed plasma 
will cross the 32 km of the geographic frame in which we have collected data in about 80 
seconds. Considerable smearing out of the plasma density is expected and an equutorward 
region of enhanced conductivity will be produced thus decreasing the calculated electric field 
there. Figure V.7 shows the effect of including this drift on the height integrated conductivi- 
ties. The microstructurc evident in the particle precipitation is now completely gone. In order 
to obtain the observed electric field, a large and fortuitously arranged divergence of the hor- 
izontal current system associated with field-aligned current densities on the order of SO /a A/m’ 
would be required. Clearly a plasma drift relative to the arc frame of this magnitude is unac- 
ceptably large for our present model. This difficulty can be overcome in two ways - an equator- 
ward arc motion or enhanced recombination coefficients. 

If the arc has a sigi.ificant eqiiatorward motion, and thi's the spcvd of the ionized com- 
ponent relative to the arc are reduced. This results in less smearing out of the ionospheric 


r 





V. ANAl.VitlS 40 

plasma density. Figure V.8 shows the height integrated conductivites resulting from various 
values for i/,, equivalent to choosing lesser values for Ey in the arc frame. Notice that as 
increases, the arc width also increases due to the fact that the data under consideration are 
obtained from the poleward going portion of the flight and any equatorward arc motion implies 
an increased arc thickness. The relative positioning of features in each case also shifts due to 
the non-linear mapping between flight time and magnetic latitude. 

Figure V.8 shows that in order to minimize the smearing effect of v,, E^ must be less than 
S mV/m. This corresponds to a of more than 300 m/s, three times the likely upper limit 
placed on the arc motion above. In addition for such rapid arc motions, the inferred arc thick- 
ness is over 100 km. Thus it does not seem that arc motion can be invoked to diminish the 
effects of plasma motion across the arc. 

If the recombination coefficients are increased we will obtain a decreased time constant for 
plasma density equilibration. This could decrease the smearing effect without the need for large 
arc velocities. A side effect, however, is that the steady state plasma density is also decreased. 
As a result the conductivities are also decreased. Equation V.l implies that for small changes 
in If, the time constant, r, for changes in the equilibrium density is approximately, 

r - — i==- (V.IO) 

2v^a 

so that to decrease t by a factor of three, « would need to be changed by an order of magni- 
tude. The equilibrium densities and conductivities would also be decreased by a factor of three. 
It does not seem likely that the height integrated conductivities should be so low (see e.g. 
Brekke ci al„ 1974), nor does it seem probable that the values used for the recombination 
coefficient are in error by an order of magnitude. Thus we are inclined to dismiss this possibil- 
ity as well. 


The conclusion reached after these considerations is that the ionized component must 
somehow have a smaller velocity with respect to the arc than would be obtained from using the 
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KxH drift. The results of our more naive model in which v, was un independent parameter are 
closer to the observed data. Perhaps the transition to the ExB drift speed occurs at a higher 
altitude than obtained from our model. This could be due to un underestimate of the neutral 
atmosphere density or the temperatures at higher altitudes, leading to an underestimate of the 
electron-neutral and ion-neutrul collision frequencies, lo the remaining discussion we wilt 
assume that, for whatever reason, v, - m, - 0 and note that u small v, (< 50 m/s) would be 
adequate to smooth the results of the simulation while any larger values would begin to elim- 
inate the observed correlation. 

We would like to now consider the elTects of v„ and on the simulation. Recall that n, 
is the y-component of the neutral wind velocity rciat.ve to the geographic frame and that it is 
used to conveniently lix the otherwise irrelevant c-component of the arc velocity. A non zeio 
value of n, is manifested by a /ero-level shift on the ob.served values of /•.*,. I'or Instance, if we 
chose i i) - 5 mVVm instead of the observed 0 mV/m in the geographic frame, a neutral 
wind speed of 


(/.\(0)-/.;<'U))) 

11 ^ « f t V . 1 1 ) 

a 

or about 100 m/s, westward, (since /\‘'(0) - 0 in this case) is implied. Recalling that neutral 
wind speeds rarely exceed 200 m/s we will limit ourselves to values of |/;\ (0)1 less than 10 
mV/m, We will also limit ourselves to values of |v„|, the v-component of the neutral wind 
speed, less than 100 m/s. This, more stringent, restriction results from our conclu.sion that the 
ioni/ed component was somehow more rigidly tied to the neutral particles and that the speed of 
the former was very small. 


I’igure V.9 displays the profile resulting from the application of our numerical pro- 
cedure to the data with ii, - v, - 0 and various values of //, and v„. The figure shows that the 
Ey profile is very scn.silivc to the choice of (0) and also fairly scivsitivc to the choice of neu- 
tral wind speed. Notice that regardless of the neutral wind speed chosen, A\(0) - -10 mV7m 



produces profiles with values of > 200 mV/m, a value rarely observed in the ionosphere. 

At this point it is necessary to ask ourselves what features constitute a good Ht to the data. 
If we settle on reproducing the observed net change in E^ of 93 mV/m (0 to -93 mV/m in the 
geographic frame) then we discover that a range of choices for (/,. and v„ is possible. The 
values are, however, dependent and the dependence may be fuirly accurately expressed by the 
relation 

Uy - lOOm/s - .6i’„ (V.12) 

within the range of acceptable values for the parameters. If, on the other hand, we concen- 
trated on reproducing the overidl f 9 rm of the observed f, profile or some portion of that form, 
a diflerent conclusion might be reached. In general «, and i’„ will be dependent on one another 
and a certain range of values for one will produce acceptable results provided the other is prop- 
erly adjusted. Since we did not tneasure the neutral wind speed we can say no more except to 
note that the observed profile can be fairly accurately reproduced with reasonable values of 
t/, and 

Finally we would like to discuss the basic source of the observed anticorrelation. As was 
pointed out in the discussion of the simple model, the inverse relationship between the electric 
field and the precipitating particles would seem to imply a linear relationship between the height 
integrated conductivities and the energy of the auroral electrons. The simulation routine was 
used to calculate the conductivities as a function of peak energy for three dilTcrent values of the 
source temperature. The source temperature is the temperature of the plasma sheet electrons 
which arc accelerated to become the auroral beam and is essentially a measure of the slope of 
the distribution function above the secondary peak. Figure V. 10 shows the results of these cal- 
culations. Notice that regardless of the source temperature, the conductivities are fairly linear 
functions of the peak energy for energies between I and 4 keV, the energies observed in this 
rather mild arc. For higher energies the Pedersen conductivity in particular levels off. This is a 
result of the increased ioni/ation at lower altitudes which docs not enhance the Pedersen 
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conductivity nearly as much as the Hall conductivity. Thereforf., for more intense arcs the 
linear behavior is non-existent and the anticorrelation would be expected to be less exact. 

V.B. On the liiah amplitude wave observations and the wave-particle effects 

In the previous section it was shown that the electron distribution function undergoes a 
rapid change after the detonation of the barium shaped charge experiment. The large secondary 
peak at downward pitch angles which was present in the pre-barium data is flattened into a pla- 
teau in the post-barium data. The introduction of the barium into the iu>.ospheric plasma 

* 

appears to have acted as some sen of catalyst for a process which removes the free energy in 
the distribution function which is represented by the secondary peak. Barium shaped charge 
experiments have been observed to stimulate lower frequency, ion cyclotron waves (Koons and 
Pongratz, 1979) due to the direct impul.se of the explosion. The observations reported here 
appear to depend only on the presence of the barium cloud and not on the explosion or plasma 
turbulence associated with the explosion. 

The question might be raised, how can the peak survive at all? Such secondary peaks 
require a region of positive slope in the distribution function which is one of the most funda- 
mental sources of plasma instability. Plasma waves which can interact with the electrons in the 
unstable region may extract energy from the particles causing a redistribution of particles from 
higher energies to lower energies and thus eliminating the peak. In the linear theory, the rate 
at which these waves convert the free energy of the particles into wave energy is proportional to 
the amplitude of the wave itself. The process is described in terms of an exponential growth 
rate, y, which depends on the wave mode and the particle distribution function. 

A particle distribution function is not necessarily unstable as a result of the existence of 
such a secondary peak. The question of instability can only be answered by detailed calcula- 
tions based on the di.spersion relations of all applicable wave modes and the entire particle dis- 
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iribution function. A great variety of deviations from a Maxwellian distribution may result in 
the growth of certain wave modes. 

In auroral arcs, the particles that these waves interact with are not a fl.(ed population. The 
auroral beam is produced continuously by the acceleration region. Thus in order for the waves 
to have an appreciable effect on the observed distribution function, they must be of sufficient 
amplitude that the power they extract from the beam is comparable to the power input to the 
particles in the acceleration region. 

In order for the waves to attain these amplitudes they must remain in contact with the 
auroral beam for some time. This requirement may be difficult to achieve since the auroral 
beam is typically confined to a latitudinally narrow volume of space (Maggs, 1976; Maggs, 
1978). For wave packets with group velocities perpendicular to the auroral strip, unless 
where w is the width of the auroral beam, the wave packets will not grow to any 
appreciable amplitude before they leave the arc and their effect on the ob.served particle distri- 
bution will be negligible. 

The restriction is not nearly as severe if the wave packets travel along the arc. Such a 
case has been treated by Maggs (1976) and the extended period of auroral contact may result in 
substantial growth for certain wave modes. The requirements on the auroral geometry however 
are such as to rule out such a mechanism ir all but a small number of real cases. 

Thus in general the beam is effectively stable by virtue of its limited spatial extent and the 
resulting brief period of contact for the waves. In the first half of the flight the particle data 
indeed reveal a stable secondary peak. The high frequency wave data indicate only occasional 
emissions near the plasma frequency and at the second electron cyclotron harmonic. Notice 
that at the times of the latter emissions the plasma frequency has dropped to a value which 
would make the upper hybrid frequency, w,?/, - + lo}, essentially the same as the second 

electron cyclotron harmonic. This fact will be important in the subsequent analysis of the 
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dispersion relation for these waves. These occasional emissions may correspond to the arrival 
in the vicinity of the ejcctable probe of wave packets which have been fortunate enough to 
remain within the auroral beam for relatively long periods. Still their existence in any one loca< 
tion is of too transient a nature to cause any noticeable redistribution of the auroral electron 
distribution. The oa'urrence of the barium detonation alters both of these observations. The 
auroral beam is rapidly eroded, eliminating the region of positive slope from the distribution 
function, and very large amplitude plasma waves are observed for an extended period. Evi- 
dently the barium cloud has somehow formed an environment conducive to energy transfer 
from the beam to the plasma waves. 




Dispersion relations 

Now we would like to examine the dispersion characteristics of waves at the observed fre- 
quencies Notice in Figure IV. I that the most prolonged emissions occur at twice the electron 
cyclotron frequency. Weaker emissions occur simultaneously at a frequency corresponding to 
the unperturbed plasma frequency. During this period 2/,,. ~ 2.4 MHz while /,,, = 2.05 MHz. 
Thus to the accuracy of our measurements the upper hybrid frequency is equal to or slightly 
less than the second electron cyclotron harmonic. This is the transition point in the first two 
branches of the dispersion relation for electron Bernstein waves, electrostatic waves traveling 
perpendicular to the magnetic field. 

Barbosa (1976) has examined the dispersion relations for high frequency electrostatic 
waves and extended the analysis to oblique wave vectors. The electrostatic portion of the 
dielectric tensor may be written 
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frequency, plasma frequency, temperature, and mass. Equation V.I3 assumes that 
|6ii| « |b» - ff I for all n. This condition is necessary to limit the effects of cyclotron damping 
on the waves of interest. Setting cn - 0, Equation V.I3 can be solved to yield 
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It is convenient to plot the dispersion relation as contour plots in a temperature norm.al- 


ized index of refraction space, d - b/u 
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Figure V.ll shows six different sets of 


constant wave frequency contours in d-space for fixed plasma frequencies in the range jf 
interest. (Note that since fl is essentially constant we consider changes in o>^ to be due to 
changes in the plasma density and hence the plasma frequency.) The solid contours in each case 
correspond to the lower frequency branch the electron Bernstein mode. The dashed con- 
tours correspond to the higher frequency branch. These plots give the contours as a function 
of positive </n and one component of di. Since the dispersion relation is an even function of </n 
the full three-dimensional iso-frequency surfaces are constructed in the obvious way by mirror- 
ing across the d, axis and rotating 360° about the d\\ axis. With this three-dimensional picture 
in mind we note the following characteristics. 

In all cases, "bulbs" form around the </|i axes which increase in size with wave frequency. 
These are just the off-axis generalization of the electrostatic plasma waves which have the 
dispersion relation 
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for ki - 0 and where w,. is tlie electron thermal speed. These bulbs collapse to the plasma fre* 
quency when kn - 0. 

In addition toroidal surfaces form around the axis when is less than twice the elec- 
tron cyclotron frequency. The size of these toroids decreases with increasing wave frequency 
and they collapse to the lesser of the upper hybrid frequency, ■■ + 1, or twice the elec- 

tron gyrofrequency, 2, when k^ - 0. 

The curves of Figure V.ll may also be used to determine the propagation direction of the 
wave packets. Since the group velocity is just the gradient of w in k-space, wave packets will 
have propagation directions given 'by the direction perpendicular to the iso-frequency surface 
toward larger frequencies. Care must be taken in this and subsequent discussion not to confuse 
motion in real space with motion in k- or d-space. 

Wave growth 

A proper calculation of the amplihcation of waves in the region of the barium cloud would 
involve several steps. Having found the dispersion relation which applies strictly only to the 
thermal background particles, we would then need to calculate the growth rates for waves of 
different frequencies as a function of the local plasma frequency and d. This would be done by 
using the measured electron distribution function including the energetic electrons which were 
considered a minor perturbation in solving for the real part of the wave frequency. A spatial 
model for the barium cloud density along with the dispersion relation and growth rates would 
then be input into a ray tracing program to compute the convective growth of waves of various 
frequencies. 

Since we have incomplete knowledge of the electron distribution function and since any 
detailed model for the barium cloud density proHle would be somewhat speculative, this pro- 
cedure does not seem warranted. Fortunately it is also largely unnecessary. Several authors 
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(Kaufman ei al„ 1978a; Kaufman et al., 1978b; Lotko and Maggs, 1979; Kaufman, 1980; Lotko 
and Maggs, 1980; and Maggs and Lotko, 1980) have shown that distributions displaying 
d//dv; > 0 are most unstable to upper hybrid waves in the lower ionosphere. Kaufman ei al. 
(1978a), Kaufman et al. (1978b), and Kaufman (1980) analyzed an observed distribution func- 
tion with a moderate secondary peak and found it weakly unstable to these waves. The pre- 
barium distribution function displayed a much more dramatic peak than the one analyzed by 
Kaufman so it seems safe to assume that an instability was present. Also as it turns out, the 
most important properties of the barium cloud are its peak density and the fact that it elongates 
along the magnetic held direction and remains fairly narrow in width. 

It is not sufficient to say that the distribution function is unstable. There will in general 
be a limited region of d-space which is unstable. Only those waves which interact more 
strongly with the electrons in the region of positive slope than with the thermal background will 
display growth. 

In Kaufman's analysis, it is shown that waves may grow when 

£ (//i „ + //,,„)> 0 (V.16a) 

where 



(V.16b) 
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(V.I6c) 
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(V.16d) 


where the argument of the Bessel functions is and vh, is the velocity of a resonant 

electron. The process of finding unstable wave modes amounts to picking a wave with a fre- 
quency bi, wave vector components, and an n such that the region of positive slope in the dis- 
tribution function lies at vh - vn, and at the value of which maximizes Jj. In addition other 
terms in the infinite sum of Equation V.16a which will generally be negative must be kept 
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small. We are generally interested in the strong emissions at the upper hybrid frequency. For 
these waves we are forced to choose n - 2. 

The main contribution to the growth of these waves comes from the perpendicular peak 
and specifically Hi 2 - In order to decrease damping from the //n terms we must keep An small. 
Looking back at Equation V.I3, recall that the approximation bn « |oi - »| was made. For bu 
violating this condition not only is cyclotron damping important, but the dispersion relation 
itself may be in error. In any case, waves with large b|, or dn are not expected to grow and we 
can limit our attention to waves close to the bi or di axis. For our purposes we will consider a 

wave to be strongly damped by the thermal electrons if d^ > ^ . Equation V.lbd also 

3(0 

limits the range of An over which we may expect growth. The unstable perpendicular peak is 
only observed at values of I’n > -6x10* cm/s and by extrapolation only exists up to a max- 
imum of I’ll ~ 2x|0’ cm/s. The most positive values of occur for |»’||| < 6x10* cm/s and 

dvi 

thus by Equation V.16d we expect growth only for values of |A,i| > LJloJ - «1 m"' where we 
have used fl - 7 5x10'’ s~' For an electron temperature of .2 eV this translates to 


It/, 1 1 > 


1 \(o - n 
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Thus the two restrictions on </|, arc compatible and yield the region 


< |,/„j < (V.17) 

3U cu 3 (i> 

in which wave growth may be expected. For du too small the waves interact with portions of 
the distribution function at large t’,i where there are few particles and no secondary peak. For 
f/ii too large, interaction with the thermal background becomes dominant. This results in strong 
damping. Another way of looking at this restriction is that large d,i implies small I'n, and a cut 
across the distribution function at small enough i,i will include the region of highly negative 

associated with the thermal electrons. 



Finally we need to consider restrictions on di . In order to maximize Hi ] we need to pick 
ki which maximizes vjO) at ~ 2x10’ cm/s ± 4x10* cm/s where the region of posi- 
tive lies. Since Jy has a maximum when the argument is 3 this restriction implies 

dvi 

ki ~ 1.1 ± .2 m"' or, 

di - .014 ± .003 (V.18) 

For large di , /j peaks at small where is highly negative. These waves are likely to be 

dvi 

damped. At small di , Ji peaks at large v’l and little or no wave growth is expected. 

» 

Wave packet trapping 

We now would like to qualitatively analyze the convection and growth or decay of specific 
wave packets as they interact with the barium produced plasma density enhancement. To do 
this we consider the dispersion relation of a given wave frequency as the plasma frequency 
(density) varies. In Figure V.12 contours of constant plasma frequency are plotted in d-space 
for a given wave frequency. Lines have been drawn indicating the regions of strong damping, 
weak damping, and growth as given by Equation V.17 and Equation V.18. 

Once again the bulb shaped contours form around the r/n axis for all values of w, but now 
they decrease in size with increasing plasma frequency collapsing to the wave frequency when 
/tii - 0. Also the toroidal surfaces are present for wave frequencies smaller than twice the elec- 
tron cyclotron frequency, but now they increase in size with increasing plasma frequency and 
when ki •" 0 they collapse to the plasma frequency that makes Z - Z„,,. Another interesting 
feature of these plots is that when 5 - 2 all the contours collapse to A i, - Ah - 0. 

Consider a wave packet of a given frequency, 5, moving through a region of space in 
which the plasma frequency. 5^, varies from to Due to refractive effects we expect 
that in time the wave will tend to visit evsry region of d-space available to it. Looking at Fig- 
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ure V.I2 we see that if S > 2 then, except tor a small range of plasma frequencies just above 
the wave frequency, the wave will have access to large values of (/« where it will be strongly 
damped. Thus we do not expect to see any waves at other than the plasma frequency above 
20. This is the situation observed on this flight. These waves are driven primarily by a secon- 
dary peak at low pitch angles. Since these particles were not measured we will not speculate 
further on the details of their occurrence except to note that the observation of plasma waves 
might be expected to coincide with pockets of diminished plasma density. An examination of 
the dispersion relation for waves near the plasma frequency shows that such pockets will be 
able to trap waves near the miniipum plasma frequency spatially and confine them to small 
values of d^. A similar mechanism will be shown to hold for upper hybrid waves except that 
the requirement in that case is that there be a small plasma density enhancement. 

If 5 < 2 there is a range of 5^ from 5^ - - 1 up to some < 2 in which f/n is 

bounded. These are the toroidal contours. Notice that when Zip - - 1 then <5 - Su/,. 

Notice also that these waves cannot exist in regions where - 1. There is a cutoff at 

this value of mp and any wave propagating toward less dense regions will be reflected at this 
point. 

Now consider a simple model of the barium cloud. It is a cigar shaped plasma density 
enhancement extended along the magnetic field lines. In the cajie of this flight the ambient 
plasma frequency was just slightly less than \!5Q. This is the plasma frequency which makes 
o)„/, - 20. The plasma frequency within the cloud rose to well above 2f) and decayed as the 
cloud expanded to the ambient value after about 40 seconds. Since the ambient plasma fre- 
quency is less than y/3ii. waves at frequencies less than 20 will have access to all portions of 
d-space within the curve corresponding to Zip^^^. When Sp^^^ is above 2 these waves will often 
take on large values of d,i and thus be strongly damped. Thus no waves should be observed 
until decays to less than 2. In fact waves are observed slightly prior to this according to 



the calculated value of Up but the rule appears to be obeyed within the uncertainty of the meas* 
urements. 

As S. continues to decrease the waves are confined to smaller and smaller values of (/« 
and the possibility for growth increases. Figure V.I3 shows detailed dispersion relations for 
w 1.97, 1.98, and 1.99. For u - 1.97, Figure V.13a shows that most waves having access to 
the region of growth are damped by virtue of the large values of </; that they have access to. 
For S 1.98, Figure V.13b shows clearly that if S. <1.74 the wave is confined to regions 
of growth or weak damping at worst. For u - 1.99, Figure V.13c shows that all waves having 
access to the growth region at this frequency are damped by virtue of the large values of ih 
which they attain at smaller values of d^. From this analysis then, we expect waves at w - 1.98 
to be the most likely candidates for convective growth. Notice that is not important for 

this consideration. In this case the curves collapse to d - 0 for 5^ - - 1 - 1.709. Even 

if is greater than this the wave is still confined to favorable values of d. To understand 

'min ^ 

the importance of ^e must consider the spatial propagation of these wave packets. If 
- 1 the wave packets are unbounded spatially and will escape the density 
enhancement and thus the auroral beam. If - 1 the wave packets are trapped 

within the density enhancement due to the reflection at the cutoff point. As a result they 
interact with the beam for a prolonged period and may achieve substantial growth. The fre- 
quencies that arc observed correspond to normal modes of the ambient plasma since it is the 
cutoff at the borders of the plasma density enhancement that traps the waves. 


Thus the most favorable conditions for growth of upper hybrid waves can be summarized 
as follows: 

1 The electron distribution function displays a highly positive at perpendicular pilch 

iiv 


angles. 
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2 The ambient plasma frequency is just under >/3in. 

3 A moderate pl(»ma density enhancement yielding Just over %/3n exists within the 
auroral beam. (Extreme enhancements tend to yield dumping due to the larger values of 
r/ii that are allowed.) 

These are precisely the conditions which were observed during this flight and thus the 
large amplitude upper hybrid waves result as a natural consequence of the fortuitous conditions. 
These waves have a wavelength of about 6 m determined by the value of at which the elec- 
tron distribution peaks. 

IVa vc-pariicle interact torn 

The observed changes in the electron distribution function are now naturally ascribed to 
the wave-particle interaction produced by the catalytic action of the cloud. Presumably at high 
altitudes (he electron distribution function retains its pre-barium peak. As the electrons enter 
the cloud where the upper hybrid wave packets are trapped, they interact with these waves. 
The waves liberate the free energy represented by the secondary peak and by the time the elec- 
trons are observed at the payload, the peak has been eroded into the plateau characteristic oi 
the post-barium distribution function. The high amplitude waves were observed for about 20 
seconds but the distribution function was changed for virtually the remainder of the flight. 
Presumably, then, wave activity continued at higher altitudes for an extended period. 

An attempt was made to calculate the rate at which energy was lost from the electrons by 
assuming the pre-barium distribution to be extant above the cic ' The results were that 
energy was conserved to the accuracy of the very rough calculations. Thus it appears that the 
waves extract little power from the beam but e.xert a large influence by redistributing the 
energy. Notice that the peak energy of the auroral beam was growing during the development 
of the barium cloud and thus such a calculation based on the assumed nature of the high alu- 



V. anai.vkis 


54 


tude distribution function is highly uncertain. 

The enhanced fluxes of low energy electrons deserves some attention at this point. Since 
the largest of these events occurs simultaneously with the highest amplitude upper hybrid 
waves, the two would seem to be associated in some way. The existence of other bursts during 
the early stages of the barium cloud and at times not correlated with wave observations may be 
further evidence that the waves are prc.ent elsewhere producing the low energy electrons which 
then travel to the observation point. These electrons have velocities of about 3x]0^ cm/s and 
will thus interact strongly with waves near the second electron cyclotron harmonic when 
ki I'l/n 3. whi'.’h once again mjiximizes This implies a of about 8 m which is about 
the value obtained by the 2.2 MHz daughter wave of the upper hybrid waves as discu.>sed 
below. Enhanced fluxes are not produced by a resonant interactioti. Electrons must be 
accelerated from a higher density, lower energy portion of phase space. This would mean that 
the interaction would have to begin at lower energies, perhaps at higher altitudes where the 
dispersion relation is ditterent. After resonantly trapping a group of low energy electrons, the 
resonance condition might shift to higher energies as the electrons move, accelerating the elec- 
trons. 

Figure V.14 is a schematic representation of the wave-particle interaction taking place 
within the barium cloud. The upper hybrid wave packets are trapped within the barium cloud 
by the cutoff which exists in the dispersion relation for just under V3. The waves modify 
the distribution function of the precipitating electrons from its form above the cloud to its form 
at the point of observation. 

M-V; ir- wa interactions 

Now we turn our attention to the 200 and 400 kHz sidebands observed near the peak in 
the growth of the upper hybrid oscillations. Due to the simultaneous appearance of oscillations 
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at the diflcrcncc frequencies these are interpreted as the product of a parametric decay process 
pumped by the high amplitude upper hybrid waves. Such a three wave interaction requires fre- 
quency and wave vector matching. That is, 

Wo ” (i>i + wi (V. 19) 

ko-k, -i-kj (V.20) 

First we will consider the requirements of frequency matching. Figure IV. 8 shows that 
Equation V.19 is satisfied for two cases, 


and 


“k/i ‘**200^ + <*»’ 2 M 

f 


(V.2Ia) 


“*W400A +Wjoi/ (V.21b) 

where the subscripts refer to the ob.servcd frequencies of th: daughter waves in Hz. The pro- 
duction of the sidebands above w,,/, is probably due to a recoupling of the low frequency waves 
to W|,/, giving the processes. 


and 


Wii/I + ‘*'2(M)A U>2 b\i 


(V.22a) 


w,,/, + W4(H)/^ — wj m/ (V.22b) 

In addition the several small oscillations that are observed at higher frequencies may also be 

generated in some manner by couplings between the upper hybrid oscillations and its high fre- 
quency daughters. The frequency of the largest of these emissions is 2w„,, - 411. This may be 
due to a coupling of the upper hybrid wave to itself. 

In order to limit the scope cf this analysis we will confine ourselves to the consideration 
of the process indicated in Equation V.21a. The next task is to find modes which exhibit the 
proper wave vector matching properties. We have no way of knowing for sure what all of the 
participating modes arc or even whether they arc electrostatic or electromagnetic. Assuming 
tnat all the modes arc electrostatic, which .seems more likely due to coupling considerations, the 
low frequency wave would be a whistler on the resonance cone and the high frequency 
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daughter would be another electron Bernstein wave wi’’^ n oblique k. The whistler has a con- 
stant k„/ki. The Bernstein wave has constant k«/ki for si. II k and becomes constant k and 
more oblique as k increases. In our specific case the 200 kHz whistler has a propagation direc- 
tion tan"'(Aii/^'i) = 12“ The 2.2 MHz Bernstein wave has a propagation direction 
tan"'(/(|i/A'i) = 55 “for k < 5m~' and eventually attains k = 9m"' for propagation directions 
within 45 “of perpendicular to B. 

As a result of the constant k,i/ki behavior of the low frequency whistler there are a multi- 
tude of possible k-matching schemes. Four possibilities are sketched in Figure V.I5. Two of 
these use the constant k portion of the Bernstein dispersion rektion and have A| = Aj »Aq. 
The other two use the constant k,i/ki portion and have Aq == A, = A 2 . In all four cases the 
pump wave is assumed to be exactly oblique. This is inconsequential in the first two cases since 
Aq is small, but in the last two cases deviations from perfect obliqueness will affect, but not 
eliminate, the matching. These four cases share the feature that B is included in the plane of 
the matching scheme. This need not be the case however it is probable that the coupling is 
stronger this way due to polarization considerations. It is also possible to allow the pump wave 
to have a propagation direction similar to the whistler in which case an electromagnetic wave 
with small A might fit the requirements for the high frequency daughter. 


^y/iy 200 kHz? 

At least one important question remains unanswered in the preceding analysis. What fac- 
tors result in the selection of 200 kHz as the preferred difference frequency? The wave vector 
matching schemes presented here could apply to any frequency up to a few hundred kHz. 

There are several possible answers to this question. First, given an instability which is 
equally likely to proceed in a number of directions, the slightest push in one of these may be all 
it takes to make the decision. Once a given decay mode is chosen it is likely to be strongly pre- 



fcrrcd due to u sell'- reinforcing stimulation. There is u radio transmitter at Luleu, Sweden, a 
distance of .^00 km from Kiruna, which operates at a frequency of 200 kik. While radiation 
from this source ut this frequency would be expected to be strongly shielded from the upper 
ionosphere, it may have played a role in determining the decay products. 

Another possibility has to do with the generation of u standing wave resonance within the 
barium cloud density enhancement. Since the waves were trapped by the cloud, it is possible 
that certain k’s might have been selected by virtue of their resonance properties. This seems 
somewhat unlikely however due to the .short wavelengths (I- 10 m) of the daughter waves com- 
pared with the si/e of the cloud (> IJ km in any dimension). 

A final interesting possibility is related to the behavior of the dispersion relation for waves 
at the upper and lower 200 kll/ sidebands of the pump. As suggested by Figure V.ll, these 
waves, having w’s of 1.83 and 2.17, when propagating in a region where <Z,, ~ y/3, have very 
similar k’s for all oblique propagation directions. Perhaps this fact, plus the common low fre- 
quency mode which couples both waves to the upper hybrid pump provides extra stimulation 
for the selection of these frequencies. 

The multiple sideband .structure may be related to this mechanism ics well. The initial 
decay into the 200 kll/ and 2.2 Mh/ daughters is reinforced by the recoupling of the 200 kll/ 
wave to the pump forming a 2.6 Mil/ daughter which shares wave vectors with the 2.2 Mil/ 
wave. The resulting 400 kll/ difference stimulates the growth of a low frequency daughter at 
this frequency which couples with a 2.0 Ml!/ Bernstein mode extracting more power from the 
pump and so on. Thus growth of waves in such a harmonic scries is self-reutlorcing due to the 
multiple couplings which arc possible in such a scheme and the process goes on until saturation 
at the point where the parametric decay matches the particle stimulated grow th of the pump. 
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VI. CONCLUSIONS 

An Aries sounding rocket was launched into a quiet auroral arc and returned comprehen- 
sive data on precipitating particles, fields, and plasma density within the arc. Two striking 
features of the data have been analyzed here. The first, the electric field-particic energy 
anticorrelation, is a large scale feature of the quiet arc. The second, the wave-particle interac- 
tion, is a small scale phenomenon artificially induced by the barium shaped charge experiment. 

The anticorrclation observed between the D.C. electric field and the peak energy of the 
auroral beam appears to be a result of the modulation of the ionospheric conductivitie • due to 
the ionization produced by the energetic precipitating electrons. The ionospheric conductivities 
increase as the peak electron energy increases due to generally increased number fluxes, larger 
ionization per incident electron, and enhanced penetration to lower altitudes where the collision 
frequency is larger. With the assumption of a negligible field-aligned divergence of the ionos- 
pheric current system, the electric field is required to vary in some sort of inverse fashion. The 
picture is complicated by the necessity to consider Hall and Pedersen conductivities, neutral 
winds, ionization flow velocities, and arc motion. A computer simulation was developed to 
properly treat the problem and the results show fair agreement with the data. As a result the 
nearly exact inverse relationship observed on this flight appears in retrospect to have been for- 
tuitous, but the inferred mechanism is more than likely operative in most quiet arcs. An appli- 
cation of the simulation routines to other rocket data obtained in quiet arcs would be helpful in 
evaluating the validity of the model. 

The wave observations appear to be due to the catalytic action of the barium plasma den- 
sity enhancement produced within the energetic electron precipitation region by the barium 
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experiment. The precipitating electron distribution function is unstable to the growth of per- 
pendicularly propagating high frequency electrostatic waves. In the unperturbed case this insta- 
bility does not lead to large wave growth because of the limited spatial extent of the arc. Waves 
move out of the arc before growing to detectable amplitudes. In the presence of a plasma den- 
sity enhancement the situation changes. Waves may exist within the enhancement that cannot 
exist outside. These waves are refractively trapped and interact with the electrons on much 
larger time scales. For the ambient plasma conditions observed on this flight, the waves 
trapped by the barium cloud were just those to which the electron distribution function was 
most unstable - the upper hybrid mode. These waves then parametrically decayed to a multi- 

I' 

tude of daughter modes producing the rich spectrum observed. As a result of the enhanced 
wave-particle interaction the free energy represented by the secondary peak was liberated and 


the electron distribution function was stabilized. 
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TABLE 1: EJECT ABLE PROBE INSTRUMENTATION 



ORIENTATION 


SENSITIVITY 


TYPK 

TO SPIN AXIS 

RANGE 

RESOLUTION 

SAMPLE TIME 

Low Lncrgy 

O'. 30*. 

6 to 500 eV 

5 7*10'* cm*stcr eV 

004 s 

Electrons 

and ISO* 

128 steps 

Af. - 30eV 

1.024 s/sweep 

High Energy 

0* 

350 eV to 30 keV 

1 8*10 *cm*sterlE) 

004 s 

Electrons 


128 steps 

A//f - 09 

1 024 s/sweep 

High Energy 

0* 

500(Z) eV to 40(/l keV 

1 1*10'* cm*stcr(f ) 

004 s 

Ions 


128 steps 

AtVt - 06 

1 024 s/.sweep 

. . 

Ion Composition 

O' 

5001/.) eV to 40(Z) keV 

1 1*10 ’ cm*sterl£'l 

.0005 s/sample 

D C Electric Field 

0“ 90° 

> 

0 to 100 mV/m 

1 mV7m 

002 s 

High Gain "D C " 

0*. w 

0 to 1 mV/m 


001 s 

Electric Field 


3 to 128 Hz 



Low Frequency 

O'. 90' 

Oto SO kHz 

300 Hz 

004 s 

Spectrum Analyzers 


512 steps 


2.048 s/sweep 

High Frequency 

0°. 90° 

0 to S MHz 

30 kHz 

.001 s 

Spectrum Analyzers 


512 steps 


.512 s/ sweep 

Magnetic Field 

0°. 90“ 

± .5 Gauss 

30 y 

004 s 

Density 


10* to lO*" enr* 


.001 s 

A C, Density 


Oto 1 kHz 

Ah/« > .01 

.0005 s 

Temperature 


500 to 4000 TC 


004 s 

A C Temjteraiurc 


Oto 100 Hz 

A r/ r > .01 

,002 s 




TABLE 2: RECORD 1 DATA 


*I>ATA 

1 UNCTION 

UNITS 

ROUTINi: 

TYi’li 

lOVI'Rf 

ln|)ui nic 11) 


SIMSIK 

Ini 

IDI/IX 

/ or profile ID 


SIMCUR 

Ini 

lUINAl 

Output lite ID 


SIMfUR 

Ini 

NCBINS 

# of cncrity bins 


SIMSK 

Ini 

IMIN 

1 owcst energy 

cV 

SIMSK 

Real 

IMAX 

Highest energy 

eV 

SIMSK’ 

Real 

TSRCi: 

Source temperature 

eV 

SIMSPC 

Real 

DSRCU 

Source density 

cm ' 

SIMSK 

Real 

HR A IK) 

Magnetic Held ratio 


SIMSK- 

Real 

nilTAl) 

Angle of B 

Degrees 

SIMI’RD 

Real 

VMl'T 

Neutral Mind s|K'ed 

m/s 

SIMDNS 

Real 

PHII) 

Angle belMcen *„ and Bo 

Degrees 

SIMDNS 

Real 

DX 

latitudinal step si/e 

km 

SIMDNS 

Real 

NXSTPS 

# of steps 


SIMSK' 

Ini 

LY 

1 ongitudinal comp of K 

mV/m 

SIMDNS 

Real 

BMACi 

B at earth's surface 

(iauss 

SIMDNS 

Real 

LI IRli; 

or /, iiidicator 


SIMl'UR 

1 og 

ICl'RNT 

/ profile indicator 


SIMCIIR 

1 og 

LIX 

, profile indicator 


SIMCUR 

Log 

lilTRJl 

Llectrojet current 

A 

SIMCUR 

Real 

vi*o\m:r 

■’article ixiwer law 


SIMSK' 

Real 







TABLE 3: RECORD 3 TO END DATA 


TmTa ’ 

FUNCTION 

’"units' 

ROUTINL 

VO 

Accelerating potential 

V 

SIMSIf 

KX 

Latitudinal comp of F 

niV7m 

SIMCUH 

CX 

Latitudinal height integrated current 

A/m 

SIMCUR 

C V 

Longitudinal height integrated current 

A/m 

SIMCUK 

I'll) 

Height integrated Pedersen conductivity 

Mho 

SIMCUK 

MALI 

Height integrated Hall conductivity 

Mho 

SIMCUK 

tlLlX 

Particle power dissipation 

erg/cm’ s 

SIMSPC 

tJOL'lt 

Field power dissipation 

erg/cni^ s 

SIMCUR 

CDRNT 

Field aligned current at 250 km 

/i A/in^ 

SIMSK' 

Q(I71) 

Production rate profile 

cm ’s ' 

SIMPRD 

IMNSd’l) 

tieciron density profile 

cm ■ ’ 

simdns 

Cl)X(171» 

Latitudinal current density profile 

M A/m^ 

SIMCUR 

CnY(l71) 

Longtiiidinal current density profile 


SIMCUR 

CD/(17I» 

Vertical current density profile 

fiA/n}^ 

SIMCUR 

SKiMAP(171» 

Pedersen couduclivity profile 

mho/m 

simcnd 

SI(iMAII(|7ll 

Hull cundudivii) profile 

mho/m 

SIMCND 

Tr<l()0» 

Downward difierentlal electron flux 

cm"^s*'eV"i 

SIMSPC 
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FIGURE CAPTIONS 








Pigure 1.1 A model of the magneioiphere. 

Figure II I A sketch of an ejectable probe and its instrumentation. (1-6) FIcctric Held and plasma booms (a) FSA 
view dircctitin (b-d) KI*A view directions, (c) 'Belly-band' antenna, (f) Magnetometer, (g) Sun sensor. 

Figure II 2 Particle deicctors. (a) An RPA. (b) an ESA 


Figure III I The trajectory of Porcupine 2 including the value of the magnetic coordinate - L. 

Figure III 2 Magnetic projection of the trajectory to 110 km 
Figure IV I Electron and plasma wave spectrog.ams 

Figure IV 2 The eguatorward, westward, and total DC electric Held The large disturbance at / - 373 seconds is due 
to the detonation of the barium shaped charge esperiment. 

Figure IV, 3 The magnetic projection to 110 km of the observed electric Held. The vectors arc labeled aa'ording to 
whether they were observed on the upward or downward leg of the llight. 


Figure IV 4 A fifty second average of the ion (luxes observed between r “ 250 and 300 seconds This .spectrum is 
typical of that observed throughout the llight 


Figure IV 5 

Figure |V 6 

Figure IV 7 

Figure IV 8, 
observations 


A comparison of the invcse total DC electric field with the peak energy of the precipitating electrons 
Sweeps of the high energy electron detector from before and after the barium detonation. 

Polar plot of the pre- and post-barium distribution function 

Twelve successive sweeps of the high fregucncy spectrum analyzer at the time of the high amplitude wave 


Figure IV 9 A detailed plot of the seventh spectrum analyzer sweep of the previous figure 

Figure V I (a and b) The relationship between the ionospheric electric field and the peak electron energy for a model 
in which the ionospheric potential configuration is calculated by adding a potential inferred from the peak electron ener- 
gy to a simple assumed magnetospheric potential configuration (c) The required magnctosphenc potential 
configuration ’.o prixluce an inverse relationship between the peak energy and the ionospheric field. 

Figure V 2 Coordinate system for the auroral simulation 

Figure V 3 The D C electric field rotated 30 "from north toward west into the coordinates of the arc 
Figure V.4 Height integrated conductivities obtained for simplest case. 
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Fifure V 5 ReUiive velociitei t)f n«uir<al wind, tin, tmi fcoirapliK frames 
Fiyure V 6 The f\ profile obtained from ibe timpfest ease compared with the observed profile 
Figure V 7 lleight imetraied conduciiviiiet In the ease u, - 0 and i, > i l9mV7m/0 
Figure V g Height integrated eonduetivities for various values of ir. 

Figure V 9 The F, profile resulting from various choices for ir, and > 

Figure V 10 Height integraicd conductivities as a function of peak energy of toe auroral beam for three different 
value* of the source temperature 

Figure V II Consiani wave frequency dispersion contours for various lived values of (a) 17 (b) s/j ici 18 Id) 
I 9 (e) 2 0 (D 2 I Curve labels are S 


Figure V 12 Constant plasma frequency dispersion contours for various lived values of S. (a) I I (b) I 9 let I 95 id) 
2.05 (e) 2 1 (0 2 2 Curve labels are Up. 

Figure V I J Details of the dispersion relation for u equal to (a) I 97 (b) 1.98 (c) I 99 

f 

Figure V 14 Diagram of barium cloud indicating proposed Interaction of trapped waves with the precipitating auroral 
electrons 


Figure V 15 
Figure VII 1 
Figure V|| 2 
Figure VII 3 
Figure VII 4 
Figure VII 5 


F'our simple wave vector matching schemes for the parametric decay process. 
A fiowchari for SIMSPC 
A fiowchari for SIMPRD 
A flowchart for SIMDNS 
A flowchart for SIMCND 

A flowchart for SIMCCR (at General (b) Detail of ihc calculitional loop 


FIgu.'c VII 6 The ditl'ereniial electron llus (TF) at three imsilions within the model arc 

Figure VII 7 Simulated arc data as a function of laiitudmal distance (a) The Input licld-aligned potential drop fb) 
The field-aligned current carried by the precipitating elecirons (c) The height inicgraied conductivities (d) The equa 
lorward component of the electric field. (e)The two comiKmenis of the height integraicd current (f)Thc particle 
power Input and Joule healing produced In ihe Ionosphere 

Figure VII 8 The projection of the i am, .• components of the current density onto the > plane 

Figure VII 9 The projection of ihe < and , components of the current density onio ihe >- j plane Downward com- 
ponems correspond to posiiive 
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Figure III.l 
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Figure III. 2 
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Figure IV, 3 
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Figure IV. 8 
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Figure V.2 
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Figure V.5 
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Figure V.14 
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